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The  stopping  power  formula  from  Bethe's  theory  contains  terms  which  are 
known  only  approximately  and  must  be  estimated  with  the  use  of  experimental 
data.  These  terms  include  a-  material  constant,  the  mean  excitation  energy  of 
the  medium,  and  the  shell-,  Bloch-  and  Barkas -corrections . In  an  analysis  of 
measured  proton  and  alpha-particle  stopping  powers  and  ranges,  modifying 
parameters  have  been  introduced  into  these  corrections,  and  the  mean 
excitation  energy  was  simultaneously  adjusted,  so  as  to  get  the  closest 
possible  agreement  with  experimental  results.  Such  an  analysis  is  reported 
here  for  elements  with  atomic  numbers  Z > 57.  The  modification  parameters 
introduced  for  the  shell  corrections  have  a simple  relation  to  atomic  energy 
levels.  The  Bethe  theory  with  the  adopted  mean  excitation  energies  and 
proposed  adjustments  of  the  shell-  and  Barkas -corrections  predicts  stopping 
powers  that  are  in  close  agreement  with  experimental  values , within  the 
experimental  uncertainties.  This  agreement  was  obtained  for  protons  with 
kinetic  energies  above  about  0.5  MeV,  and  for  heavier  ions  of  charge  z at 
energies  above  (z  - 1.5)  MeV/u. 


Report  of  work  done  under  contract  for  the  National  Institute  of  Standards  and 
Technology 
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I.  INTRODUCTION 


The  Bethe-Bloch  theory  of  stopping  power,  S,  (Bethe, 
1930;  Bloch,  1933;  Fano,  1963)  with  a Barkas  effect  correction 
term  and  including  shell  corrections  has  been  shown  to  agree 
well  with  experimental  data  for  protons  and  alpha  particles 
with  kinetic  energies  T/M  (M  is  the  number  of  nucleons  in  the 
particle)  above  0.5  MeV  traversing  light  elements  for  which  I- 
values  had  been  determined  independently.  A single  free 
parameter  for  the  Barkas  effect  was  introduced  (Bichsel  and 
Porter,  1982).  For  heavier  elements,  semiempirical  I-values 
and  shell  corrections  were  used  to  calculate  S for  comparison 
with  experimental  data  (e.g.  Bichsel,  1961;  Bichsel,  1972; 
Janni,  1982;  Porter  and  Bryan,  1984;  ICRU,  1984).  Since  shell 
corrections  for  the  M-shells  (Bichsel,  1983)  and  the  separate 
L“subshells  (Bichsel,  1987)  are  now  available,  I have  explored 
the  validity  of  the  theory  with  these  functions.  It  was 
necessary  to  introduce  corrections  for  the  outer  (N,  0,  P,  Q) 
shells:  they  were  derived  from  the  M-shell  corrections  with  a 
scaling  procedure  using  atomic  structure  data.  For  all 
elements  with  57<Z<92,  the  I-value  was  the  only  free  parameter 
not  simply  related  to  Z.  The  theoretical  Bloch  function  and 
an  empirical  function  for  the  Barkas  effect  (Bichsel,  1990) 
were  used. 

The  present  study  is  an  outgrowth  of  earlier  work 
(Bichsel,  1961,  1964,  1967)  and  is  an  effort  complementary  to 
those  by  Andersen  and  Ziegler  (1977),  Ziegler  (1977)  and  Janni 
(1982).  Many  observations  about  problems  with  experimental 
data  were  made  in  these  references  and  should  be  studied 
there . 


II.  THEORY 

The  theory  of  S for  heavy  elements  is  complex  (see  Table 
XI  below).  It  is  not  possible  to  calculate  S a priori  with  an 
accuracy  of,  say,  1%,  and  empirical  modifications  of  the 
current  a priori  theoretical  functions  are  needed  to  get 
calculated  values  agreeing  with  experimental  data  S^. 
These  modifications  are  attained  by  varying  the  values  of 
parameters  in  the  functions.  Many  of  the  parameters  are  quite 
interdependent:  a change  in  a parameter  causing  an  increase 
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in  S over  an  extended  energy  range  can  readily  be  compensated 
by  a change  in  one  of  the  other  parameters  causing  a decrease 
in  S (see  Table  IV  below). 


It  is  the  purpose  of  the  present  paper  to  consider  the 
complete  theory  of  S outlined  below,  and  parameters  used 
should  have  values  plausibly  related  to  atomic  data.  In  other  » 
analyses,  incomplete  theories  (e.g  Janni,  1982,  did  not 
explicitly  use  z and  z corrections)  or  polynomial  fits  to 
the  shell  corrections  (where  the  parameters  have  no  physical 
significance)  were  used  (Andersen  and  Ziegler,  1977).  The 
same  goal  is  achieved  with  all  these  approaches: 
semiempirical  functions  are  given  which  approximate  the 
experimental  data. 


Further  improvements  in  the  theory  of  interactions  of 
charged  particles  with  matter  have  been  achieved  for  collision 
cross  sections  (e.g.  Anholt,  1979,  Cohen  and  Harrigan,  1985), 
but  similar  improvements  have  not  been  made  for  the  stopping 
power.  Relativistic  corrections  for  atomic  properties  are 
discussed  in  the  Appendix,  but  were  not  used  in  the  data 
analysis 


A.  Stopping  power,  S. 


The  expression  used  for  the  calculation  of  the  stopping 
power  of  fast  charged  particles  is 


S 


^ _ JL  Z 2 
dx  ” ^2  A ^ 


[1] 


with  T the  kinetic  energy  (MeV)  of  the  particle, 

X the  absorber  thickness,  g/cm'^, 

k = 47re^N  /mc^  = 47rr^  N me  = 0.307072  MeV  cm^ , 
o'  o o 

V - (3c  the  speed  of  the  incident  particle, 

7^  = 7=  1 + T/M^c^, 

ze  the  charge  of  the  incident  particle, 

e the  electron  charge, 

m the  electron  rest-mass,  mc^=51 1,000  eV, 
c the  speed  of  light, 

r = 2.817941*10”  cm  the  classical  electron  radius, 
o 23  / 

N = 6.022045*10  atoms/mole,  Avogadro ' s number, 
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Z the  atomic  number  of  the  absorber, 

A the  atomic  weight  of  the  absorber  (in  g), 

L the  stopping  number  (B  in  older  papers),  and 
M the  rest  mass  of  the  particle;  p:  M^c^  * 938.2561  MeV 
° a:  MqC^  = 3727.316  MeV 

Either  /?  or  T will  be  used  as  the  variable  indicating  the 
particle  energy,  with  /?^« ( T/M^c^ ) • ( 2+T/MqC^ ) / ( 1 +T/MqC^ ) ^ . 

For  particles  heavier  than  electrons,  the  stopping  number 
L is  expressed  in  the  form 

L(z)  s + zL^  + ^2^^^ 
with 

L (/?)  = t(0)  - j2n  I - + (Giz.P)  - 6{0))/2  [2a] 

where  ti/3)  ~ J2n  (2mc^  7^ ) - /?^,  I is  the  mean  excitation 

energy  of  the  absorber,  C(/?)  the  total  shell  correction, 
G(z,/?)  the  Mott  correction  term,  6 the  correction  for  the 
density  effect,  L.|  the  Barkas  correction  term  and  the  Bloch 
correction  term.  It  is  useful  to  define  the  experimental 
value  of  the  stopping  number,  L^,  obtained  by  solving  Eq.  ( 1 ) 
for  an  experimental  value  of  the  stopping  power: 

L^(/?)  = / (kz^z).  [2b] 


B»  I“ValueSo 

The  mean  excitation  energy  I is  defined  by 

00  00 

fin  I = J f(E,0)  fin  E dE  / J £(E,0)  dE,  [3] 

0 0 

00 

J f (E,0)  dE  = 1 , 

0 

where  E is  the  energy  transfer  in  a transition  and  f(E,0)  the 
dipole  oscillator  strength  DOS  (Fano,  1963),  which  is  related 
to  the  optical  absorption  coefficient  (e.g.  Barkyoumb  and 
Smith,  1990).  The  Bloch  parameter  is  defined  by 
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b = I/Z  [3a] 

For  the  Thomas-Fermi  atom,  Bloch  ( 1 933  ) showed  that  b is  a 
constant . 

For  many  gases,  I-values  have  been  determined  with  Eq. 
(3)  (see,  e.g.,  Zeiss  et  al.,  1977;  Jhanwar  et  al.,  1983). 
The  only  metals  for  which  I was  calculated  with  Eq.  (3)  are 
Al  (Shiles  et  al.  1980)  and,  partially,  Si  (Bichsel,  1 988  ). 
For  heavy  elements,  insufficient  information  is  available 
about  f(E,0)  to  permit  a determination  of  I. 

For  metals,  the  excitation  function  for  the  collective 

excitation  of  the  valence  electrons  should  be  used  in  Eq.  (3) 
( "plasmon  excitations",  e.g.  Raether,  1980).  Qualitatively, 
the  influence  of  plasmons  may  be  understood  from  the  value  of 
the  most  probable  energy  loss  hf^  for  valence  electrons,  given 
in  Table  VI;  hf^  was  derived  from  electron  energy-loss  spectra 
given  by  Ahn  et  al.  (1  983  ).  hf^  depends  on  the  structure  of 
the  metal  and  the  number  of  valence  electrons  and  thus  will 
not  relate  simply  to  Z. 

A quantity  related  to  I is  the  plasma  energy  hf  (Fano, 
1963),  defined  by  (hf)^  = 830.4  p Z/A  (hf  in  eV,  p the  density 
of  the  metal  in  g/cm^ ) , shown  in  Table  I.  This  quantity 
replaces  the  I-value  in  S for  very  high  particle  energies. 

Clearly,  hf  depends  strongly  on  the  density,  explaining  the 

small  values  for  Pb,  Bi  and  Th.  The  small  values  of  hf^  and 
hf  for  Pb  compared  to  those  for  Au  may  explain  in  part  why  the 
experimental  value  of  the  Bloch  parameter,  b^,  for  Pb  is 
smaller  than  that  for  Au. 

Another  approach  to  getting  information  about  stopping 
power,  I-values  and  shell  corrections  is  to  use  the 
statistical  atomic  model  (Bonderup,  1967).  Values  of  I for 
atoms  were  calculated  by  Bichsel  and  Laulainen  (1971)  with 
this  model  using  the  relativistic  wavef unctions  of  Liberman  et 
al.  (1965,  1971).  The  resulting  Bloch  parameters,  b^,  based 

on  7=1.347,  are  given  in  Table  I.  Similar  results  were 
obtained  by  Chu  and  Powers  (1972).  Values  for  the  solids 

would  differ  by  several  percent;  for  light  elements  there  are 
large  differences:  for  Al,  b(atom)  = 9.56  eV,  Dehmer  et  al. 

(1975);  b( solid ) = 1 2 . 77  eV,  Shiles  et  al.  (1980);  for  Si, 
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b(  atom ) *9 . 39  eV,  b(  solid  ) =*  1 2 . 43  eV,  Bichsel  (1  988).  It  must 
also  be  clearly  understood  that  the  DOS  used  in  this  model 
differs  strongly  from  actual  optical  data  (Johnson  and 
Inokuti , 1983). 

It  is  tempting  to  try  to  derive  I-values  from  the  well 
known  experimental  ionization  energies  (e.g.  Bearden  and 
Burr,  1 967  ).  I have  used  an  expression  given  by  Sternheimer 
et,al.  (1984)  ((their  Eq.  (8));  it  was  derived  for  a 
different  purpose)  to  calculate  I from  J^,  the  plasma 
frequency  hf  and  with  an  estimate  of  DOS  for  photon  energies 
above  the  ionization  energies  expressed  by  a factor  r . Values 
bg  calculated  with  a constant  value  r=2  are  given  in  Table  I. 
These  results  demonstrate  that  one  should  resist  this 
particular  temptation. 

The  experimental  results  for  b^  from  the  present  study 
are  given  in  Table  I,  and  are  discussed  in  section  III  B. 
Experimental  values  used  by  Andersen  and  Ziegler  (1977),  Janni 
(1982)  and  ICRU  (1984)  are  also  given. 

C«  Shell  corrections  C(v,z). 

7.  General  review. 


In  the  Bethe  theory,  shell  corrections  must  be  calculated 
on  a shell-by-shell  basis: 

C(v,Z)  = ECi,(v,Z)  = Cj^(v,Z)  + Cj^(v,Z)  + Cj^(v,Z)  + ...  [4] 

V 


Usually,  the  dependence  on  particle  speed  v is  expressed  in 
terms  of  the  variable  r? : 


= (mv^/2)  / mv^-Z^/2j  = (mv^/2)  / 


[4a] 


with  e^=RZ^,  v^  the  Bohr  speed  .(v^=c/137),  R=mv^/2  = 13.6  eV 
the  Rydberg  energy  of  the  hydrogen  atom,  and  Z^  the  effective 
charge  of  the  absorber  atoms  for  electrons  in  shell  The 
dependence  of  on  the  atomic  number  Z also  enters  via  the 
ionization  energy  J^,  expressed  in  terms  of: 
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[4b] 


'^V  ' ! «!/• 

In  principle  (Janni,  1982),  these  functions  should  be 
calculated  for  each  subshell  in  the  atom  (Bearden  and  Burr, 
1967).  Here,  a somewhat  simpler  approach  was  used,  as 

outlined  below. 

Shell  corrections  for  K-  and  L-shells  were  derived  by 
Walske  (1952,  1956)  with  the  nonrelativistic  hydrogenic 
approximation.  Recently,  the  corrections  for  the  M-shells 
(Bichsel,  1983)  and  the  L-subshells  have  been  calculated 
(Bichsel,  1987),  also  with  the  hydrogenic  approximation. 

For  the  outer  shells  no  calculations  have  been  made,  and 
a scaling  procedure  was  used  in  which  it  was  assumed  that  the 
shell-corrections  for  outer  shells  have  the  same  shape  as 

those  for  the  inner  shells  ( Hirschf elder  and  Magee,  1948; 
Bichsel,  1961;  Janni,  1982).  Vertical,  V^,  and  horizontal, 
H^,  scaling  factors  were  introduced: 

where  v stands  for  any  one  of  the  outer  shells  and  /x  stands 
for  the  inner  shell.  Presumably,  would  be  proportional  to 
the  ratio  of  ionization  energies,  and  V related  to  the 

number  of  electrons  in  the  shell.  In  order  to  assess  the 
plausibility  of  this  approach,  it  is  instructive  to  compare 
the  shapes  of  the  known  shell  corrections.  In  Fig.  1,  they 

are  shown  for  K,  L and  M electrons  in  gold.  The  functions 

have  been  plotted  in  such  a way  as  to  coincide  at  the  maximum 
value.  Clearly,  the  functions  are  similar  in  shape,  and  the 
scaling  procedure  can  be  used  with  some  justification,  but  the 
horizontal  scaling  factors  f^  given  in  the  Fig.  are  only 
approximately  proportional  to  the  ionization  energies  W^,  and 
the  vertical  scaling  factors  g^  are  not  well  correlated  with 
the  number  n of  electrons.  At  present  it  is  not  known  how 
closely  the  hydrogenic  calculations  approximate  the  correct 
functions.  Because  empirical  parameters  H and  V are  used  in 
Eq.  (5),  the  degree  of  approximation  cannot  be  assessed  from 
experimental  data  for  S.  More  accurate  calculations  are 
needed. 
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Bonderup  (1967)  derived  shell  corrections  from  the 
stopping  power  of  an  electron  gas  (Lindhard  and  Winther,  1964) 
with  electron  density  varying  with  distance  from  the  nucleus 
(’’local  plasma  model”)  and  gave  functions  based  on  the  Lenz- 
Jensen  model.  Empirical  functions  for  the  shell  corrections 
were  given  by  Andersen  and  Ziegler  (1977)  and  by  Janni 
(1982).  Shell  corrections  for  Au  from  these  sources  are  shown 
in  Fig.  2. 


2.  Present  approximation . 

(a) ^  Inner  shells  (K^  L,  M) 

The  functions  given  by  Walske  (1952),  Khandelwal  (1968), 
and  Bichsel  (1983,  1987)  were  used.  The  effective  charge, 

in  Eq.  (4a),  was  assumed  to  be  nf  , where  n is  the  principal 
quantum  number  and  f the  orbital  exponent  given  by  dementi  et 
al.  ( 1 967  ).  For  the  K-shell,  Zj^=Z-0.3  was  used.  The 

quantities  and  for  some  elements  are  given  in  Table  II. 
From  the  variations  in  for  the  subshells  and  from  the 

differences  in  shape  seen  in  Fig.  1 it  is  evidently  advisable 
to  calculate  for  each  inner  subshell  separately. 

(b) .  Outer  shells. 

From  the  values  of  the  scaling  factors  g^  and  f^  in  Fig. 
1 we  must  conclude  that  for  the  outer  shells  there  are  no 

compelling  choices  for  the  scaling  factors  and  in  Eq. 

( 5 ) . They  must  be  chosen  such  as  to  give  good  agreement  with 
experimental  data. 

The  ionization  energies  for  the  outer  shells  of  some  of 
the  heavy  elements  are  given  in  Table  III.  In  order  to  reduce 
the  number  of  free  parameters,  and  in  view  of  the  uncertainty 
of  the  choice  of  it  is  reasonable  to  consider  shell 

correction  functions  for  groups  of  subshells  for  which  the 

ionization  energies  are  similar:  to  Njjj/  and  N^,  and 

all  the  other  shells.  I have  used  this  approach  and  have 

named  the  functions  ^N3 ' six  parameters  , 

, H2,  and  . For  the  first  and  second  group,  the 

number  of  electrons  in  the  shells  to  is  constant  (18), 
and  we  would  assume  that  the  factors  and  should  change 
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slowly  with  Z,  at  most.  It  appears  plausible  to  use  initial 
values  of  given  by 


H 


u 


[5a] 


where  n^  is  the  number  of  electrons  in  subshell  u (Table 
III).  I found  that  values  and  calculated  with  this  Eq. 
gave  better  agreement  with  experimental  data  for  Au  than 
values  different  by  ±20%,  and  no  further  searches  were  made 
with  different  and  Thus  there  were  now  only  the 

parameters  and  for  which  values  were  undetermined.  It 
was  expected  though  that  and  should  be  near  1 . A further 
option  was  the  choice  of  the  inner  shell  m from  which  the 
outer  shell  functions  are  scaled.  I have  found  the  best 
agreement  with  experiments  by  using  the  functions  for 

'^N2'  ^MV  Cjj,  and  0^3- 


For  all  electrons  outside  of  (in  gold  a total  of  33 
electrons,  with  ionization  energies  between  0 and  108  eV ) a 
single  function  obtained  by  scaling  with 

and  as  free  parameters  (H^  might  be  or  greater,  and 

might  be  expected  to  depend  on  Z ) . 


The  method  of  determining  the  parameters  and 

is  described  in  section  IIH.  It  will  be  seen  that  these 
four  parameters  are  ample  to  provide  calculated  values  of  the 
stopping  power  agreeing  well  with  experimental  data. 

By  choosing  the  shell  corrections  to  fit  experimental 
data,  they  will  also  compensate  for  errors  in  and  z . 

D.  Barkas  and  Bloch  corrections,  and 
3 4 

The  need  for  z and  z corrections  in  the  theory  was 
established  experimentally  by  Andersen  et  al  (1977).  Basbas 
(1984)  discussed  the  problem  in  a general  context.  Bichsel 
( 1990  ) analysed  experimental  data  and  found  that  for  Au  only 
an  empirical  function  for  approximated  experimental  data 

(extending  from  1 to  4 MeV ) well: 
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[6] 


(/?)  = 0.002833  •/?  ' 

Here,  this  expression  has  been  used  for  all  proton  energies, 
and  for  all  Z>57.  These  extrapolations  are  dubious. 

The  term  L2  is  written  in  the  form  derived  by  Bloch 
(1933): 

00 

Lj  = '('(■>)-  Re[i;'(  1+iy )]  = -y^  ^ 2~ FJ 

j = l 

where  ^ is  the  logarithmic  derivative  of  the  T function,  and 
y * 2v  /v  a 'icxlP  (a=1/l  37.036  is  the  fine  structure 
constant).  For  y <1,  the  sum  is  equal  to  1.202...  does 

not  depend  on  Z,  but  does  include  terms  of  all  even  powers  of 
Zo  Possible  errors  in  L.^  and  will  be  compensated  by  the 

choice  of  the  shell  correction  parameters. 

E.  Charge  state  corrections  and  nuclear  collisions. 

Bichsel  (1990)  found  that  it  was  not  necessary  to  use  any 
charge  state  corrections  for  protons  and  ot  with  T/M>0.5  MeV  (M 
is  the  number  of  nucleons  in  the  ion).  For  Li-ions,  a reduced 
charge  z S3  appears  to  be  needed  for  T/M<2  MeV,  but  no 
definite  form  of  2 could  be  derived  from  the  experiments. 
For  lower  energies  and  heavier  ions,  the  need  for  charge  state 
corrections  will  be  seen  in  the  Figs.  and  in  Table  VIII. 

Nuclear  collisions  contribute  less  than  0.1%  to  the 
stopping  power  at  the  energies  considered  here,  and  thus  are 
neglected. 


F.  Mott  correction  G and  density  effect, 

Ahlen  (1978,  1980)  gave  a close-collision  correction 
G(z,/?)  due  to  the  Mott  cross  section.  For  p and  a,  a rough 
approximation  is  G « z/?/137,  where  -10<a^<2  for  0.3<T/M<30 
MeV.  Thus  this  term  is  less  than  0.1%  for  p and  a,  but  will 
be  more  important  for  heavy  ions  (e.g.  about  2%  for  3000  MeV 
Ca  ions ) . 
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The  algorithm  used  here  for  the  density  effect  has  been 
given  by  Sternheimer  et  al.  (1  984).  It  must  be  noted  that 
approximations  were  made  which  may  introduce  errors  of  the 
order  of  0.3%  into  at  T>100  MeV  (see  Figs.  6.2  to  6.4  in 
ICRU-37,  1 984  ).  The  effect  amounts  to  less  than  0.04%  for  20 
MeV  protons  in  gold,  and  thus  is  only  important  for  some  of 
the  measurements  at  high  energies  and  for  the  heavy  ions. 

G.  Ranges  and  multiple  scattering. 

Ranges  were  calculated  from  stopping  power  S with  the 
csda  ( continuous-slowing-down ) approximations 

T 

R(T)  = + J / S(T')  [8] 


For  present  purposes,  for  protons,  T^  = 0.4  MeV,  and  R(T^)  was 
taken  as  the  total  pathlength  given  by  Janni  (1982).  For  10 
MeV  protons  in  Au,  this  contribution  amounts  to  less  than  1% 
of  the  total  range.  Due  to  multiple  scattering,  experimental 
projected  ranges  are  shorter  than  R(T).  Corrections  for  this 
effect  were  given  by  Bichsel  and  Uehling  (1960),  Bichsel 
(1960),  Berger  and  Seltzer  (1964),  Bichsel  (1972),  Janni 
(1982),  and  by  Bichsel  and  Hiraoka  (1989). 

H.  Method  for  parameter  determination. 

The  parameters  to  be  determined  are  the  I-value  and  the 
scaling  factors  and  V^  of  Eq.  (5)  for  the  outer  shell 
corrections.  They  were  found  from  experimental  data  with  a 
least-squares-deviation  procedure . 

In  earlier  studies  (e.g.  Porter  and  Bryan,  1984),  the  I- 
value  was  introduced  explicitly  as  a free  parameter  in  the 
data  analysis.  This  can  be  avoided  with  the  following 
approach.  Customarily,  the  deviation  6 is 
calculated,  and  the  sum  H 6^  is  considered  as  a function  of 
the  five  parameters.  A best  fit  is  obtained  if  the  sum  is  a 
minimum.  A variation  of  the  parameter  I can  be  avoided  if  the 
equation  for  Si/S)  is  rearranged  as  follows: 


[9] 


6{0)  + Inl  = Y(/?)  = t^(0)  - c^(0)  - + z L^(|3) 


+ 1.^(0)  - L^{0) 

where  t^{0)  = t(0)  + (G(z,/?)  - S(0))  / 2 


c^(0) 


Co(/?) 


C^{0)  + C^(^)  + C^(/?) 


and 


the  dependence  on  Z is  not  indicated » This  expression 

implicitly  contains,  in  c^(/?),  the  four  free  parameters 
of  Eq.  (5).  If  there  are  no  systematic  errors  in  the  theory, 
we  expect  that  Y(/?)  will  be  a constant  (but  subject  to 
stochastic  errors  of  the  data),  and  expiYiP) is  an 
experimental  I-value  for  each  data  point.  We  define  as  the 
average  of  p experimental  values  of  Y(/?): 

Y-  = E / p [10] 

^ p 

and  by  assuming  Y =lnl,  the  average  deviation  a defined  by 

a 

= E / 'P"”  • 

is  the  same  as  2 ^ / (p-1)-  With  this  approach,  the 

parameter  search  is  performed . for  a space  reduced  by  one 
dimension  (i.e.  a five  parameter  search  is  reduced  to  a four 
parameter  search).  depends  on  the  parameters  H^,  and 

can  be  considered  as  a function  of  the  f our-dimensional  space 
with  coordinates  V^.  The  smallest  define  sets  of 

parameters  giving  best  fits  to  experimental  data. 

2 

It  was  found  that  there  were  many  local  minima  of  cr  , and 
therefore  the  method  of  steepest  descent  was  not  suitable  for 
the  parameter  determination;  a grid  search  was  used  instead. 
The  parameters  for  small  were  recorded,  and  the  associated 
functions  Y(/?)  were  plotted  versus  particle  speed  /?  and 
examined  for  systematic  deviations.  If  the  deviations 

YiP))  were  randomly  distributed,  or  if  they  deviated 

systematically  by  much  less  than  the  experimental  error, 
satisfactory  values  of  the  four  parameters  and  the 
corresponding  I-value 


1 1 


exp(Yj 


[10b] 


I 


a 


had  been  found.  I is  subject  to  systematic  errors  of  S . 

a X 

III.  EXPEEIMENTAL  DATA. 

For  most  elements,  only  the  experimental  data  for  protons 
were  used  for  the  parameter  determination.  The  data  for  ot- 
particles  in  Au  were  used  implicitly  by  the  determination  of 
the  function  for  the  Barkas  effect  (Bichsel,  1990).  For 

measurements  below  30  MeV,  relative  to  A1  or  Cu  (e.g.  Burkig 
and  MacKenzie,  1957),  theoretical  values  from  Bichsel  (1972) 
were  used  to  calculate  S^. 

I have  found  that  systematic  differences  occur  between 
experimental  data  from  different  sources  (see  the  Figs.)  and 
that  the  uncertainties  assigned  by  the  authors  to  their 

experimental  data  represent  only  a qualitative  measure  of 
This  means  that  a simultaneous  analysis  of  several  data 

sets  is  not  practical.  Therefore,  best  fit  parameters  were 
determined  for  each  data  set,  and  average  values  of  the 
parameters  were  then  used.  For  some  data  sets  where 

stochastic  errors  were  less  than  1%  (e.g.  Semrad,  1990; 
Oberlin  et  al.,  1980,  1982;  Matteson  et  al.,  1978),  a smooth 

function  obtained  from  a three  parameter  fit  (H^,  and  I) 

has  been  used  to  represent  the  data  in  the  Figs. 


For  gold,  I designed  an  average  experimental  data  set  for 
proton  energies  above  0.3  MeV.  Between  0.3  and  1.5  MeV,  data 
by  Luomajtrvi  (1979),  Semrad  (1990),  Andersen  and  Nielsen 
(1981)  and  Santry  and  Werner  (1981)  were  used  and  given  about 
equal  weight.  Between  1.5  and  3 MeV,  the  data  by  Andersen  and 
Nielsen  (1981)  were  reduced  by  0.8%.  Above  3 MeV,  the  data  by 
the  Nara  group  (Ishiwari  et  al.,  1 984;  Shiomi  et  al.,  1 986  ) 
and  those  by  S^^rensen  and  Andersen  ( 1 973  ) were  weighted 
inversely  with  their  quoted  errors.  Other  data  were  not 
used. 


For  each  experimental,  S^(T),  and  theoretical, 
value,  the  relative  difference  r(T)  was  determined: 

r(T)  - (S^^T)  - ! ^x^^^ 


St(T), 

[11] 


1 2 


Values  of  r(T),  in  percent,  are  plotted  in  the  Figs.  The 
average  standard  deviation  for  an  experimental  data  set 

with  p values  was  defined  by 

= Y.  r(T)  ^ / (p-1  ).  [12] 

r(T)  and  c should  be  compared  with  the  errors  of  given 
by  the  autiiors,  shown  in  Tables  VI,  VIII,  and  in  the  Figures. 
If,  in  general,  |r(T)|  is  less  than  ^ consider  the 

theory  to  be  adequate.  If  |r(T)|  exceeds  1^^  I there  may  be 
systematic  errors  of  the  theory  or  the  experiments,  and  the 
reader  is  invited  to  choose  which  to  believe.  Symbols  used  in 
plotting  the  data  are  given  in  Table  V. 

IV e RESULTS  OF  PARAMETER  SEARCHES. 

A.  Shell  correction  parameters  and  I-values. 

In  principle,  the  parameters  could  be  different  for  each 
element.  After  a preliminary  four  parameter  search  for 
several  elements,  I found  that  V^  and  V^  could  be  assumed  to 
be  constant  for  all  Z.  This  assumption  may  have  to  be  changed 
as  more  accurate  data  become  available.  From  an  examination 
of  the  experimental  data  I concluded  that  only  for  gold  there 
were  enough  data  to  permit  a meaningful  four  parameter 
search.  For  other  elements,  a two  parameter  search  for  and 
V^  was  made,  and  it  was  possible  to  find  fixed  values  of 
which  were  valid  for  groups  of  elements.  Only  for  V^  was  it 
necessary  to  assume  a dependence  on  Z.  Eventually,  the  I-value 
resulting  from  Eq.  (10b)  was  the  only  truly  free  parameter. 

With  the  search  program  based  on  the  use  of  Eq.  (9), 
uncertainties  cannot  readily  be  stated  for  the  parameters  (or 
the  stopping  power  function  ) . It  will  be  seen  from  Figs. 
6-13  that  this  is  not  a major  problem. 

The  four  parameters  V^ , V^ , and  V^  for  the  outer  shell 
corrections,  Eq.  (5),  were  determined  for  the  average  data 
set  for  gold  with  the  search  procedure  outlined  in  section  II 
H.  The  grid  search  for  local  minima  of  , Eq.  (10a),  was 
performed  with  20  to  30  values  of  each  of  the  four  parameters, 
with  the  initial  value  of  about  I.S-H^,  the  final  value 


about  the  initial  equal  to  0.5  and  the  final  values 

up  to  3.  In  each  search,  up  to  500,000  *p  values  of  were 

calculated.  Results  of  the  search  with  and  from  Table 
III  are  given  in  Table  IV.  The  4 parameters  V,,  v^,  V3  and  H3 
giving  a local  minimum  of  are  listed,  together  with  the 

I-value  and  (Eq.  (12)).  Since  for  all  these  parameter 

sets  the  average  standard  deviation  <7^  between  experimental 
and  calculated  values  is  less  than  even  the  smallest  quoted 
experimental  error  (±0.3%),  they  all  could  be  used  to 

approximate  the  average  data  set  well.  Evidently,  this  data 
set  is  insufficient  to  determine  the  parameters  uniquely  (this 
of  course  is  also  true  for  the  individual  sets).  Because  the 
number  of  electrons  associated  with  is  larger  than  that  for 
, a value  larger  than  is  desirable.  Thus,  the 

parameters  V^=1.25,  V2=1.4  are  henceforth  used  for  all 
elements,  and  the  value  should  be  valid  for  elements 

neighbouring  gold  in  the  periodic  table. 

It  is  expected  that  the  parameters  and  associated 

with  electrons  in  the  outermost  shells  to  Q)  will  depend 

on  the  atomic  number  Z.  Therefore  I made  a two  parameter 

search  for  and  for  best  fits  for  all  elements.  I found 

that  for  Z>73  a value  H2*13  gave  satisfactory  fits;  of  course, 

the  values  and  I varied  for  each  experimental  data  set. 

Values  of  V^  are  plotted  in  Fig.  3:  a tendency  toward  an 

increase  of  V^  with  Z can  be  discerned,  especially  for  the 

data  from  Denmark  and  Nara.  Since  V^  should  be  related  to  the 

number  of  outermost  electrons,  n =Z“46,  the  function 

o 

V^=(Z”46)/25  [13] 

was  chosen  to  represent  this  parameter.  The  value  of  V^  is 
quite  sensitive  to  experimental  uncertainties  because  at  1 

MeV  contributes  only  about  10%  to  the  stopping  number  L,  Eq. 
(2),  and  only  1 % at  6 MeV. 

For  the  elements  with  57<Z<73,  the  results  of  a two 

parameter  search  for  and  V^  are  shown  in  Table  VI.  Note 

that  <7^  for  some  data  exceeds  o slightly.  This  may  mean  that 

the  authors  underestimate  o fsee  the  Figs.).  No  definite 

e 

trend  of  with  Z can  be  seen,  but  it  appears  that  a larger 
value  is  appropriate  for  Z<60  than  for  Z>60,  thus  H2=50  and 


}i^^25  were  chosen  respectively.  For  Z=73,  a smaller  value  is 
indicated,  and  the  value  for  Au  (H^=13)  was  used.  Then,  a one 
parameter  search  was  made  for  . Values  giving  best  fits  are 
shown  in  Fig.  4 and  in  Table  Via,  together  with  and  cr^. 
Note  that  may  change  little  while  may  change  much.  This 
means  that  the  change  in  and  V^  is  compensated  partly  by 
the  change  in  In  Fig.  4,  no  systematic  trend  of  V^  with 

Z can  be  seen  below  or  above  Z=60,  therefore  I chose  constant 
values,  approximately  equal  to  the  average  for  all  data,  viz. 
V2=*3.85  for  Z<60,  ^^=^2,3  for  Z>60. 

Finally,  the  average  value  I for  each  data  set  was 
determined  with  the  parameters  defined  above  and  the  values 
b *I  /Z  are  shown  in  Fig.  5.  Again  the  <7^  are  less  than  cr 
with^ these  I^  for  most  data  sets  (Table  Via),  and  in  generaf 
they  are  only  slightly  larger  than  c for  the  one  and  two 
parameter  fits.  The  fluctuations  in  I^  for  each  Z evidently 
are  in  part  an  expression  of  the  systematic  differences 
between  experimental  data,  but  because  the  values  were 
obtained  from  data  at  different  energy  ranges,  they  may  also 
indicate  possible  problems  in  the  assumptions  about  the 
theory. 


B.  Average  I-values  for  the  elements. 


For  the  elements  with  a single  data  set  (Sm,  U ) , the 

analysis  is  finished:  the  experimental  value  of  I for  the 

element  is  I =I  . For  elements  with  several  sets  of 

e a 

experimental  data,  we  must  select  a value  I^,  a weighted 
average  of  I^  (I  by  the  definition  of  Eq.  (2a)  is  a property 
of  the  material).  Greater  weight  was  assigned  to  values  from 
higher  energies.  The  values  selected  are  shown  in  Table  I. 
They  are  only  valid  in  the  context  of  the  other  parameters 

selected  here,  and  of  the  experimental  data  used  in  the 

analysis.  An  uncertainty  of  I^  of  ±1.5%  to  ±5%  should  be 
assumed.  A comparison  with  other  experimental  values  of  I is 
only  meaningful  to  show  trends  with  Z. 


It  is  notable  that  (except  for  W,  Bi  and  Th ) b and  the 
value  calculated  with  the  local  density  model  differ  by 

less  than  ±5%  (the  average  deviation  is  (0.6±3.3)%),  even 
though  b^  was  calculated  for  single  atoms,  while  b^  was 


measured  for  the  metal.  Some  of  the  variation  in  b must  be 

e 

related  to  the  variation  in  hf  and  in  the  plasmon  energy  hf^ 
of  the  valence  electrons.  In  particular,  this  may  explain  the 
difference  in  b^  between  Au,  Pb  and  Bi.  It  appears  advisable 
to  be  suspicious  of  the  values  of  b^  for  W,  Bi  and  Th. 
Fluctuations  in  b for  neighbouring  Z may  be  indicative  of 
systematic  errors  in  the  measurements. 

The  differences  between  b^  and  bj  and  b^  are  explained  in 
part  by  the  fact  that  a larger  set  of  data  was  used  here.  It 
must  be  understood  clearly  that  values  b^  will  change  as 
further  experimental  data  become  available. 

The  I-values  for  Pb  and  U for  the  proton  energies  below 
30  MeV  are  less  than  those  indicated  from  the  higher 
energies.  Further  studies  appear  to  be  needed. 

A determination  of  I-values  with  an  uncertainty  of  less 
than  2%  from  other  methods  (e.g.  Shiles  et  al.,  1980  ) is 
highly  desirable. 

V.  Comparison  of  theory  and  experiment. 

A.  Protons. 

The  comparison  between  experimental  data  and 
theoretical  functions  for  T<30  MeV  is  made  in  Figs.  6-13. 
The  relative  difference  r(T),  Eq.  (11),  is  shown  as  a 
function  of  kinetic  energy  T/M  of  the  protons.  The  authors* 
experimental  uncertainties  are  shown  at  only  a few  energies. 
There  is  no  evidence  to  invalidate  the  theory.  It  is  somewhat 
surprising  that  the  theory  agrees  with  experiments  at  energies 
T well  below  1 MeV.  There  are  no  general  trends  for  r(T)  to 
be  definitely  larger  than  0,  thus  no  need  for  a reduced 
effective  charge  z is  evident.  For  Ta,  Fig.  7,  the 
agreement  between  theory  and  experiments  as  well  as  between 
different  experiments  is  poor.  While  the  data  by  Luomajflrvi 
(1979)  for  p in  Ta  differ  considerably  from  theory,  they  agree 
quite  well  for  W and  Au. 

Some  experimental  data  not  included  in  the  Figs.  are 
considered  next.  For  Th,  data  by  Teasdale  (1949),  Sonett  and 
MacKenzie  (1955)  and  Burkig  and  MacKenzie  (1957),  with  H2=13, 
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V^s*1.76  result  in  1^*766  eV,  <7^=±0.7%.  Energy  loss 
measurements  for  protons  with  initial  energies  between  1 4 
and  25  MeV  in  thick  absorbers  were  made  by  Bichsel  and 
TschalSr  in  1966  with  the  method  described  by  TschalSr  (1967) 
and  Tschaiar  and  Bichsel  ( 1 968 ) . The  residual  energies 
were  measured  with  silicon  detectors  and  corrected  for 
multiple  scattering.  A theoretical  thickness  t^^RC )-R( ) 
was  obtained  from  range  energy  tables  calculated  with  Eq. 
(8).  The  uncertainty  of  the  energy-measurements  (less  than 
±0.1%  for  T^,  about  ±24  keV  for  T^ ) corresponds  to  an 
uncertainty  of  no  more  than  0.5%  in  calculated  thicknesses 
t^.  The  results  have  not  been  published  so  far  because  the 
theory  used  in  1967  was  inadequate.  Experiment  and  present 
theory  are  compared  in  Table  VII.  For  both  metals,  a 
systematic  trend  with  energy  is  seen.  Larger  I-values  would 
provide  somewhat  better  agreement. 

Sakamoto  et  al.  (1989)  measured  the  stopping  power  for 
73  MeV  protons  in  ten  elements  with  an  uncertainty  of  ±0.7%. 
Data  are  shown  in  Table  VIII.  Only  the  value  for  Pb  exceeds 
(7^  significantly.  A value  1=855  eV  would  be  needed  to  achieve 
r=0  for  Pb.  For  Al,  r=0.3%  and  for  Cu,  r=0.8%  from  the  tables 
described  below. 

Several  measurements  with  high  energy  protons  have  been 
reported.  They  are  listed  in  Table  IX.  Usually  the  range  of 
particles  with  a fixed  energy  was  measured  in  Al  or  Cu.  Then 
a given  layer  in  the  principal  absorber  was  replaced  by 
another  material,  and  the  range  measurement  was  repeated.  The 
ratio  of  thicknesses  of  materials  causing  the  same  energy  loss 
was  thus  determined.  In  order  to  obtain  the  stopping  power 
for  the  second  absorber,  that  of  the  first  absorber  must  be 
known.  Therefore,  stopping  power  tables  for  Al  and  Cu  were 
calculated  with  the  shell  corrections  given  in  ICRU-37  (1984), 
with  I(A1)=166  eV,  I(Cu)=322  eV  and  the  density  effect  of 
Section  II  F.  No  corrections  for  straggling,  multiple  and 
nuclear  scattering  and  nuclear  reactions  were  made  beyond 
those  made  by  the  authors.  It  must  be  noted  that  the  change 
in  the  multiple  scattering  and  the  straggling  in  the  second 
absorber  will  lead  to  changes  in  the  total  range  (of  the  order 
of  0.5%)  which  have  not  been  considered  so  far  (Bichsel  and 
Hiraoka,  1989). 


It  was  found  that  the  ratio  of  absorber  thicknesses  even 
for  fairly  large  energy  losses  (e.g.  T^a340  to  T^  = 270  MeV ) 

differs  by  no  more  than  0.1%  from  the  ratio  of  stopping  powers 

at  the  mean  energy  (i.e.  300  MeV).  This  can  readily  be 

confirmed  from  existing  range  tables  (Bichsel,  1972). 

Therefore,  was  calculated  from  the  ratio  of  thicknesses 

("mass  stopping  power")  given  by  the  authors  and  the 

theoretical  stopping  power  of  the  reference  material  (A1  or 
Cu ) at  the  mean  energy  T given  by  the  authors. 

For  Bakker  and  Segre  (1951),  Tj^=340  MeV,  T^P^70  MeV, 
<7^*±1%,  was  calculated  with  of  both  A1  and  Cu,  and  the 
average  is  given  in  the  table.  For  U,  exceeds  c 

slightly.  For  the  ratio  S^(  A1 ) /S^(  Cu ) , the  di^erence  of  it 
may  be  indicative  of  systematic  errors  in  theory  or 

measurement . 

Barkas  and  von  Friesen  (1961)  with  750  MeV  protons 
measured  thicknesses  of  several  materials  relative  to  the 
equivalent  thicknesses  of  Cu  at  several  depths  in  a copper 
absorber.  Theoretical  values  of  the  stopping  power  of  Cu  were 
used  to  obtain  for  Pb,  U and  Al.  Since  the  estimated 

experimental  error  is  considerably  less  than  1%,  the  relative 
differences  r between  theory  and  experiment  are  disturbingly 
large.  Larger  I-values  would  be  needed  to  get  r=0.  If  the 
measurements  had  been  made  relative  to  Al,  the  average  r for 
Pb  would  be  0.5%,  for  U,  0.3%. 

Vasilevskii  and  Prokoshkin  ( 1 967  ) measured  stopping  power 
of  Pb  relative  to  Cu  for  p,  d and  The  data  are  listed  in 
Table  IX  at  the  equivalent  proton  energy.  The  absolute  values 
of  given  for  Cu  by  these  authors  (shown  in  Table  IX)  were 
used  to  calculate  for  Pb.  For  Pb,  the  agreement  between 
theory  and  experiment  is  good,  but  the  experimental  values  for 
Cu  exceed  the  theoretical  values  by  0.8%  on  the  average. 

Vasilevskii  et  al.  (1969)  measured  the  stopping  power  of 
Pb  relative  to  that  of  Al  with  the  method  used  by  Barkas  and 
von  Friesen  (1961)  with  an  initial  beam  energy  of  660  MeV. 
The  theoretical  stopping  power  function  for  Al  was  used  to 
calculate  for  Pb.  Here,  on  the  average,  for  Pb  exceeds 
S by  1.3%.  It  is  interesting  to  include  the  experimental 
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data  for  Cu:  in  contrast  to  the  data  by  Barkas  and  von 
Friesen,  there  is  essentially  no  difference  between  theory  and 
experiment . 

From  the  data  at  energies  above  20  MeV,  a higher  I-value 
(which  would  reduce  ) is  indicated  for  Pb  and  U.  The  need 
for  further  correction  terms  in  the  theory  must  also  be 
considered. 


B.  Helium  and  heavier  ions. 

Selected  data  for  a are  shown  in  the  Figs.  and  in  Table 
VIII.  The  values  of  r generally  increase  with  decreasing 
T/M<0.5  MeV,  indicating  the  need  for  a charge  state 
correction . 

Since  the  Barkas  correction  was  determined  from  the 
data  for  Au  by  Andersen  et  al.  (1  977),  we  can  expect  the 
agreement  between  S^  and  S^  for  a to  be  good.  This  is  indeed 
the  case.  Agreement  within  experimental  errors  was  also  found 
for  Li-ions  (Bichsel,  1990). 

For  Au,  only  the  data  by  Fontell  and  Luomaj^rvi  ( 1 979  ) 
(dashed-dotted  line)  and  by  Matteson  et  al.  (1978)  (solid 
line  near  T/M=s0.5  MeV)  are  shown  in  Fig.  6a.  Other  data  are 
shown  in  Table  VIII. 

For  the  experimental  values  by  Anthony  and  Lanford  (1982) 
for  C-ions  in  gold,  a charge  state  correction  is  needed  if  it 
is  assumed  that  is  correct  for  z>3o 

The  experimental  data  by  Datz  et  al.  ( 1 977  ) appear  to 
have  large  systematic  errors.  The  need  for  a charge  state 
correction  or  for  changes  in  etc.  appears  for  z>4. 

Experimental  data  by  Ishiwari  et  al.  ( 1 971  , 1 977,  1 978) 
and  data  by  Takahashi  et  al.  (1983)  generally  agree  with 
theory.  For  the  C-ions,  and  L2  amount  to  8%  and  5%  of  L, 
respectively.  The  expression  given  by  Eq.  (6)  seems  to  be 
valid  for  z = 6 at  this  energy.  The  Mott  term  amounts  to  less 
than  0.2%.  For  Pb,  an  I-value  of  747  eV  would  give  better 
agreement  with  experiments. 
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The  data  for  ^^0,  and  ^^Ca  ions  traversing  Ta 
and  Au  agree  well  with  calculated  values,  except  for  0-ions 
with  T<100  MeV  in  Au,  where  a charge  state  correction  appears 
to  be  necessary.  The  data  at  75  MeV  seem  to  be  wrong. 

A function  correcting  for  charge  state  z*  is  needed  for 
T/M<(z-1.5)  MeV.  The  data  presented  here  are  not  sufficient 
to  arrive  at  a quantitative  description  of  z {0). 

C • Kange  measurements . 

Range  data  were  not  suitable  for  the  determination  of  the 
shell  correction  parameters.  In  principle,  I-values  could  be 
obtained  from  them.  In  Table  X,  experimental  ranges,  R^,  are 
compared  with  calculated  ranges  R^,  Eq.  (8). 

The  ranges  of  protons  in  Au  measured  by  Bichsel  et  al. 
(1957),  were  corrected  by  the  multiple  scattering  corrections 
given  by  Bichsel  and  Uehling  ( 1 960  ).  They  are  in  good 
agreement  with  the  theory.  Asymmetries  in  the  range- 
straggling function  (Lewis,  1952;  Tschallr,  1968;  Bichsel  and 
Hiraoka,  1989)  have  not  been  taken  into  account  yet. 

The  ranges  for  protons  in  lead  measured  by  Bloembergen 
and  van  Heerden  (1951)  i(T  ±0.5%  to  ±1%)  are  given 
uncorrected  for  multiple  scattering.  They  exceed  the 
calculated  values  on  the  average  by  (0.2±0.3)%.  Since  the 
multiple  scattering  correction  amounts  to  ?«1.7%  (Bichsel, 
1972),  the  calculated  values  are  too  small  by  this  amount.  An 
I-value  of  860  eV  would  be  needed  to  give  calculated  ranges 
agreeing  with  measured  values,  corrected  for  multiple 
scattering.  It  must  be  noted  that  the  multiple  scattering 
corrections  given  by  Janni  (1982)  are  about  50%  larger. 

The  ranges  of  protons  in  lead,  measured  by  Mather  and 
Segr^  (1951),  corrected  for  multiple  scattering  (Bichsel, 
1960;  Berger  and  Seltzer,  1964)  agree  well  with  calculated 
values.  The  ranges  measured  in  Al  and  Cu,  though,  are  1%  and 
1.5%  less  than  the  theoretical  ones. 

A range  measurement  by  Vasilevskii  and  Prokoshkin  (1967) 
for  620  MeV  protons  in  Pb  relative  to  the  range  in  Al  exceeds 
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the  calculated  range  by  1.5%. 

The  corrected  ranges  for  750  MeV  protons  in  Pb  and  U 
measured  relative  to  the  range  in  Cu  by  Barkas  and  von  Friesen 
(1961)  exceed  the  calculated  ranges  by  1.2%.  The  range  for 
A1  is  also  given.  It  agrees  well  with  R^,  in  contrast  to  the 
data  for  the  stopping  powers  in  Table  IX. 

We  see  that  for  three  out  of  four  sets  of  range 
measurements  for  Pb  a larger  I-value  is  indicated  than  from 
the  stopping  power  data. 

VI.  USE  OF  THE  THEORY  FOR  LOW  ENERGIES. 

The  expression  for  the  stopping  number  L of  Eq.  (2) 
contains  several  terms  which  change  quite  rapidly  with 
particle  energy  at  small  speeds.  This  can  be  seen  in  Table 
XI,  where  the  terms  of  Eqs.  (2)  and  (2a)  are  shown  for 
protons  in  gold.  The  shell  corrections  are  combined  into 
"inner  shells",  c^  and  "outer  shells"  c^  (Eq.  (9)).  The  high 
speed  approximation  is  defined  by  Lg  * f(/?)  “ i2n  I,  usually 

called  the  "Bethe  approximation".  For  E>1  MeV,  L_  differs  by 
no  more  than  10%  from  L,  but  the  sum  of  the  corrections  still 
amounts  to  1.5%  at  100  MeV  (the  Mott  term  G and  the  density 
correction  S amount  to  less  than  0.1%  and  are  included  in 
f(^)). 

With  decreasing  energy,  the  various  correction  terms 
begin  to  contribute  increasing  amounts  to  L.  Around  1 MeV,  the 
net  contribution  from  the  shell  corrections,  c^+c  , is  almost 
zero.  Below  about  0.5  MeV,  the  major  contribution  to  L is 
from  c.!  Clearly,  the  values  for  L below  about  3 MeV  depend 
strongly  on  the  values  of  the  parameters  used,  and  it  is  quite 
remarkable  that  most  experimental  data  agree  well  with 
theoretical  values  at  small  energies  for  protons  as  well  as 
for  a. 

All  the  terms  in  Eqs.  (2)  and  (2a)  have  a well  defined 
physical  meaning  even  at  the  smallest  energies  listed  here, 
and  therefore,  no  definite  energy  can  be  given  at  which  the 
theory  is  invalid.  On  the  other  hand,  the  magnitude  of  each 
of  the  terms  is  not  well  determined  theoretically.  If  for 
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example  different  functions  were  used  for  and  changes 
in  and  could  compensate  for  most  changes.  For  heavier 
ions,  deviations  between  theory  and  experiment  would  be 
accommodated  by  a charge  state  correction  which  would  require 
further  free  parameters. 

VII.  TABLES  OF  STOPPING  POWER. 

Values  of  stopping  power  for  p and  ot  calculated  with  the 
present  theory  are  given  in  Tables  XII  and  XIII.  They  are 
compared  with  other  tabulations  in  Fig.  14.  Andersen  and 
Ziegler  (1977)  published  an  evaluation  of  experimental 
stopping  power  data  for  protons.  For  energies  greater  than  1 
MeV,  they  used  six  free  parameters  to  calculate  the  stopping 
power  (see  their  Table  1,  p.16:  the  parameter  A^-S . 099  • 1 Z 
is  not  a free  parameter;  A^  = 2mc^  / I is  a free  parameter  and 
corresponds  to  I^  in  the  present  paper;  the  other  five 
parameters,  Ag  to  show  a systematic  dependence  on  Z,  but 
are  not  related  to  the  parameters  and  used  here).  For 
E<1  MeV,  four  parameters  were  used. 

Janni  (1982)  gave  stopping  powers  for  protons.  He  used 
the  scaling  procedure  of  Eq.  (5)  to  obtain  shell  corrections 
for  each  subshell  from  Walske’s  (1956)  L-shell  corrections. 
Data  from  the  tables  by  Williamson  et  al.  (1  966)  are  also 
shown  in  Fig.  14.  Data  for  p in  Pb  above  1 MeV  in  Bichsel 
( 1 972  ) differ  by  less  than  1%  from  present  values.  A major 
reason  for  the  relatively  large  differences  below  2 MeV  is  the 
inclusion  in  the  present  study  of  the  data  by  Knudsen  et  al. 
(1980),  Sirotinin  et  al.  (1984)  and  LuomajSrvi  (1979)  (see 
Table  VI)  which  were  not  available  for  the  earlier 
evaluations . 

For  all  the  tables,  the  differences  seen  in  Fig.  14  are 
in  part  due  to  differences  in  the  choice  of  the  I-values  (see 
Tables  I and  X),  but  there  are  also  considerable  differences 
in  the  shell  corrections  (see  Fig.  2). 

A comparison  of  the  proton  stopping  power  tables  of 
Andersen  and  Ziegler  (1977)  and  the  tables  for  ot  particles  by 
Ziegler  (1977)  was  made  by  Bichsel  (1990). 
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VIII.  CONCLUSIONS. 


Initially,  the  theory  presented  here  contained  seven  free 
parameters:  and  z/a1,2,3  and  I , in  addition  to  the 

choice  of  for  the  L-  and  M-shells  (see  Table  II).  After 
preliminary  searches,  and  were  chosen  according  to  Table 
III,  while  and  were  determined  in  parameter  searches 

from  the  averaged  data  for  protons  traversing  gold,  and  were 
assumed  to  be  constant  for  the  other  elements.  Finally  it  was 
found  that  constant  values  of  for  groups  of  elements,  with 

given  as  a function  of  Z,  and  I^  as  the  only  free  parameter 
gave  good  agreement  with  experiments  for  protons,  a-particles 
(Figs.  9-13)  and  some  heavier  ions  (Table  VIII)  for  all 

elements  with  Z>57.  For  gold,  I estimate  the  uncertainty  of 
the  theory  to  be  ±1%  below  3 MeV,  ±0.5%  between  3 and  20  MeV, 
1%  above  20  MeV.  For  other  elements,  it  may  be  ±2  to  3%  below 
3 MeV,  ±1%  above.  The  influence  of  the  uncertainty  of  the  I- 
value  must  be  added.  For  T>30  MeV,  further  correction  terms 
may  be  needed. 

The  present  approach  is  plausible  insofar  as  it  includes 
all  of  the  elements  of  current  thoughts  about  Bethe-Bloch 
theory.  It  relies  heavily  on  the  experimental  data  for  the 
determination  of  the  parameters  though.  In  examining  the  data 
for  individual  experiments  for  various  Z in  Tables  VI  and  VIA 
it  is  seen  that  a single  experimental  data  set  cannot  be 
expected  to  provide  the  parameters  H^,  V^  and  I^  suitable  for 
other  data.  Thus  it  appears  inadvisable  to  determine 
parameters  on  the  basis  of  data  for  a single  data  set, 
especially  if  they  extend  over  a restricted  energy  range;  and 
earlier  studies  (e.g.  Bichsel,  1967;  Porter  and  Bryan,  1980) 
are  only  applicable  for  the  particular  sets  of  experimental 
data  used  in  their  analysis. 

If  the  theory  is  to  be  used  for  elements  not  listed  in 
Table  I,  an  I-value  must  be  chosen.  This  could  be  done,  e.g., 
by  using  b from  a neighbouring  element  or  a calculated  value 
of  b^. 

Many  systematic  errors  of  unknown  magnitude  are 
associated  with  the  theoretical  functions  used  here.  Examples 
are: 
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nonrelativistic  hydrogenic 


wave- functions 


1 . use  of 
( Appendix ) , 

2.  use  of  scaling  procedure  for  the  calculation  of  the  shell 
corrections , 

3.  extrapolation  of  empirical  to  lower  energies  and 

different  Z, 

4.  influence  of  the  approximations  used  by  Bloch  in  his 
derivation  of 

5.  neglect  of  higher  terms  in  the  Born  approximation, 

6.  approximations  used  for  the  Mott  term  and  the  density 
effect, 

7.  neglect  of  charge  exchange  effects  (e.g.  Arnau  et  al., 

1990). 


In  addition,  systematic  errors  of  the  experimental  data 
cannot  necessarily  be  discerned-  An  example  is  the 
modification  of  the  Andersen  et  al.  ( 1 977  ) data  for  Au  by 
Andersen  and  Nielsen  (1981).  Therefore  the  results  for  I 
based  on  a single  set  of  data  (i.e.  Ce,  Pr,  Sm,  Ho,  Th,  U? 
must  be  considered  to  be  tentative.  Furthermore,  the  results 
for  Ta  (Fig.  7)  do  not  inspire  much  trust  in  the  experiments 
or  the  theory.  Only  the  results  for  Au  and  maybe  Pb  can  be 
considered  more  than  tentative.  I would  be  surprised  though 
if  new  measurements  would  show  the-  need  for  changes  in  the 
basic  parameters  of  Tables  II  and  III,  and  in  and  . 

Measurements  for  proton  energies  between  0.5  and  6 MeV  for 
several  elements  with  an  uncertainty  of  no  more  than  0.3% 
would  demonstrate  the  Z-dependence  of  and  more  clearly. 
Similar  measurements  would  be  required  to  establish  values  of 
H^,  and  I for  compounds. 

Further  developments  of  the  theory  appear  to  be  tedious 
and  may  not  be  worthwhile  unless  further  measurements  show  a 
need  for  them.  A better  approximation  for  could  be 

determined  by  measurements  at  0.5<T/M<2  MeV  at  least  for 

protons  and  a-particles  with  an  error  of  less  than  0.3%. 

The  independent  determination  of  I-values  from  Eq.  (3) 
with  an  uncertainty  of  less  than  ±2%  is  desirable  (it  would 

help  e.g.  with  the  problems  with  the  data  for  Ta).  Accurate 

X-ray  absorption  measurements  and  electron  energy-loss  data 
similar  to  those  by  Ahn  et  al.  (1  983  ) would  be  needed  for 
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this  purpose.  Such  I-values  would  also  help  in  establishing 
the  errors  of  the  scaling  factors  and  - or  they  might 
demolish  the  present  approach. 

While  approximations  better  than  the  nonrelativistic 
hydrogenic  calculations  for  K-  and  L-shell  excitations  have 
been  made  for  collision  cross  sections  (reviewed  e.g.  by 
Cohen  and  Harrigan,  1985),  these  calculations  still  differ  by 
large  amounts  from  experimental  data  at  low  particle  speeds. 
It  is  thus  an  open  question  whether  corresponding  calculations 
for  the  shell  corrections  would  be  helpful. 

It  should  be  explored  whether  the  differences  between 
theory  and  experiments  for  stopping  powers  and  ranges  at 
energies  above  30  MeV  indicate  the  need  for  further  correction 
terms  in  the  theory. 
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APPENDIX:  RELATIVISTIC  EFFECTS  IN  ABSORBER  ATOM. 


Leung  (1989)  described  a correction  to  the  Bethe  theory 
which  he  obtained  from  considering  relativistic  effects  for 
the  atomic  electrons.  For  Au,  he  estimated  an  increase  of  S 
of  about  2%  due  to  this  effect.  The  change  was  related  to 
relativistic  corrections  to  the  Bethe  sum  rule  described  by 
Rustgi  et  al.  (1988).  The  introduction  of  this  correction 
term  into  the  present  analysis  [with  the  expression  given  in 
Leung's  Eqs.  (13)  and  (14)]  changes  the  coefficients  H and  V 
for  outer  shell  electrons,  and  decreases  the  I-values.  For  an 
average  data  set  for  Pb,  consisting  of  the  data  by  Ishiwari  et 
al.  ( 1984)  and  S^^rensen  and  Andersen  ( 1973),  the  following 
three  parameter  best  fit  was  obtained:  H2  = 12,  1*720 
eV,  with  cr^s±o.08%.  The  fit  is  as  good  as  that  shown  in  Fig. 
7.  If  these  parameter  values  are  used  for  all  the  high  energy 
data  for  Pb,  the  average  value  of  r is  (0c2±0.9)%,  compared  to 
(1.4±0.9)%  for  the  standard  parameters.  Thus  the  Leung 
correction  to  the  stopping  power  function  brings  the  theory 
into  closer  agreement  with  experiment  at  the  higher  energies, 
and  it  appears  desirable  to  explore  this  effect  in  more 
accurate  studies.  Also,  more  accurate  measurements  at  both 
low  and  high  energies  would  be  useful  to  assess  the  accuracy 
of  this  correction. 
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TABLE  I.  Values  of  Bloch  parameters  b*I/Z  from  various 
sources,  and  related  quantities  (section  IIB):  average 
experimental  plasmon  energies,  hf^  and  plasma  energies  hf,  all 
in  eV.  was  calculated  for  atoms  with  a local  density 
model,  b^  is  the  experimental  value  from  Table  V and  Fig.  5 ; 
and  bg  was  calculated  from  ionization  energies.  Values  used 
in  other  tabulations  are:  bj  (Janni,  1982);  b^  (Andersen  and 
Ziegler,  1977);  b^  (ICRU,  1984).  They  differ  from  b^  partly 
because  different  shell  corrections  and  different  experimental 
data  sets  were  used,  partly  because  different  assumptions 
about  the  dependence  on  Z were  used. 


z 

hf 

V 

hf 

^L 

^e 

^A 

57 

La 

45.9 

8.53 

8.32 

11.8 

9.75 

8.42 

8.8 

58 

Ce 

- 

48.2 

8.6 

8.76 

- 

9.75 

8.5 

8.8 

59 

Pr 

48.5 

8.67 

8.64 

11.2 

9.76 

8 .59 

9. 1 

62 

Sm 

50.6 

9.04 

9.05 

10.0 

9.78 

8.84 

9.78 

64 

Gd 

51  .5 

9.08 

8.83 

9.0 

9.27 

8.81 

9.2 

66 

Dy 

53.7 

9.4 

9.17 

9. 1 

9.8 

9.09 

9.8 

67 

Ho 

54.5 

9.36 

9.55 

- 

9.8 

9.3 

9.8 

68 

Er 

- 

55.6 

9.6 

9.56 

8.6 

9.11 

9.41 

9.7 

70 

Yb 

“ 

48.3 

9.65 

9.66 

7.6 

9.82 

9.46 

9.82 

72 

Hf 

- 

66.2 

9.76 

9.32 

“ 

9.83 

9.47 

9.83 

73 

Ta 

- 

74.6 

9.78 

10.05 

9.2 

10.11 

9.37 

9.8 

74 

W 

25 

80.3 

9.8 

10.53 

9.6 

10.17 

9.36 

9.8 

77 

Ir 

86.5 

9.87 

10.23 

” 

9.93 

9.55 

9.8 

78 

Pt 

23 

84.2 

9.96 

10.08 

9.8 

10.58 

9.73 

10.1 

79 

Au 

24 

80.2 

1 0 

1 0 

9.9 

10.21 

9.56 

1 0 

82 

Pb 

15 

61  . 1 

9.79 

9.5 

10.7 

9.99 

9.26 

1 0 

83 

Bi 

18 

56.9 

9.71 

00 

0 

00 

- 

9.87 

9.22 

9.9 

90 

Th 

- 

61  .3 

9.06 

8.51 

- 

8.18 

9.17 

9.4 

92 

U 

77.4 

9.16 

9.09 

14 

9.56 

9.21 

9.7 
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TABLE  II.  Parameters  used  for  the  calculation  of  shell 

corrections  for  the  L-  and  M-shells  for  some  heavy  elements. 

n^  is  the  number  of  electrons  in  each  subshell.  The  total 

number  of  electrons  included  in  these  shells  is  26.  In  the 

first  line  for  each  element,  the  value  of  the  atomic  unit  for 

2 

each  subshell,  defined  by  Eq.  (4a),  is  given  in 

keV.  In  the  second  line,  the  ionization  energy  is  given  in 
atomic  units,  where  is  the  measured 

ionization  energy  found  in  Bearden  and  Burr  (1967). 


shell 

4 

4l 

^III 

^I 

^II 

M 

III 

M 

IV 

My 

2 

2 

4 

2 

2 

4 

4 

6 

element 

La 

23.94 

37.91 

37.91 

19.57 

19.78 

19.78 

25.22 

25.22 

57 

0.262 

0.155 

0.145 

0.07 

0.061 

0.057 

0 .034 

0.033 

Gd 

30.32 

48.48 

48.48 

25 . 98 

26.51 

26.51 

34.38 

34.38 

64 

0.276 

0.164 

0.149 

0.072 

0.064 

0.058 

0.035 

0.034 

Ta 

39.55 

63.97 

63.97 

35.5 

36.49 

36.49 

48.05 

48 .05 

73 

0.295 

0.174 

0.154 

0.076 

0.068 

0.06 

0 .037 

0 .036 

Au 

46.34 

75.51 

75.51 

42.29 

43.73 

43.73 

58 .36 

58.36 

79 

0.31 

0.182 

0.158 

0.081 

0.072 

0.063 

0.039 

0.038 

Pb 

49.66 

81  .63 

81.63 

45.99 

47.58 

47.58 

63.9 

63.9 

82 

0.319 

0.186 

0.16 

0.084 

0.075 

0.064 

0.04 

0.039 

U 

60.21 

99.65 

99.65 

56 

58.49 

58.49 

80.4 

80.4 

92 

0.361 

0.21 

0.172 

0.096 

0.085 

0.071 

0.045 

0.044 
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TABLE  III.  Ionization  energies  (Bearden  and  Burr,  1961) 
and  horizontal  scaling  factors  and  for  the  calculation, 
Eq.  (5),  of  the  shell  corrections  (the  average  value  for 
the  8 electrons  in  shells  N_  to  N__-. ) and  (for  the  10 
electrons  in  shells  and  N^).  Values  of  and  are 
defined  by  the  average  value  of  the  ratios  weighted 
with  number  n^  of  electrons  in  each  subshell. 


shell 

My 

Ni 

^II 

^III 

^IV 

«1 

«2 

6 

2 

2 

4 

4 

6 

element 

La 

832 

270 

J 

t/ 

206 

(eV) 

191 

99 

99 

3.96 

8.4 

Gd 

1 185 

376 

289 

271 

141 

141 

4 

8.4 

Er 

1 409 

449 

366 

320 

1 77 

1 68 

3.94 

8.22 

Ta 

1735 

566 

465 

404 

241 

229 

3.84 

7.43 

Au 

2206 

759 

644 

545 

352 

334 

3.6 

6.47 

Pb 

2484 

894 

764 

644 

435 

413 

3.43 

5.89 

U 

3552 

1 441 

1273 

1 045 

780 

738 

3.01 

4.71 

36 


TABLE  IV.  Parameters  V , and  V«  for  local  best  fits  to 

the  average  data  set  for  protons  in  gold,  obtained  in  a four 
parameter  grid  search.  H =3.6  and  H-=6.47  are  given  in  Table 
II.  is  defined  in  Eq.  (10a),  in  Eq.  (12),  and  I = 

exp(Y  ) [Eq.2(10)].  A local  best  fix  is  defined  be  a lol:al 

minimSm  of  in  the  four  dimensional  space  defined  by  Vw  V 
V-  and  . Note  that  the  smallest  values  of  do  not 

necessarily  ^ccur  at  the  same  grid-point  as  the  ^smallest 
values  of  , This  is  because  of  the  nonlinear  relation 
between  Hn  I and  S. 


^2 

«3 

^3 

<7^ 

la/eV 

a % 

X 

1 c3 

1 . 1 

12.8 

1 .35 

0.01494 

792 

0,1525 

1 .3 

1.15 

13. 

1 .35 

0.01466 

791  .4 

0,1555 

1 . 3 

1 .2 

1 3 . 

1 .35 

0.0151  1 

790.3 

0. 1602 

1 c 3 

1 .25 

13.2 

1 .35 

0.01631 

789.7 

0. 1686 

1 . 3 

1 .3 

13.4 

1 .3 

0.01645 

790.9 

0,1678 

1 .25 

1 .35 

12.8 

1 .3 

0.01498 

791  ,3 

0.1633 

1 .25 

1 .4 

13. 

1 .3 

0.01429 

790.7 

0.1571 

1 .25 

1 .45 

13,2 

1 .3 

0.01459 

790 

0.1613 

1 .25 

1 .5 

13,4 

1 .3 

0.01579 

789.3 

0 . 1 737 

1 .2 

1 .55 

12.8 

1 .25 

0.01674 

791  .7 

0, 1 684 

1 .2 

1 .6 

13. 

1 .25 

0.01513 

791 

0.159 

1 .2 

1 .65 

13.2 

1 .25 

0.01448 

790.3 

0,1588 

1 .2 

1 .7 

13.2 

1 .25 

0.01447 

789.3 

0 . 1 677 

1 .2 

1 .75 

13.4 

1 .25 

0.01532 

788.6 

0.172 

1 .2 

1 .8 

13.6 

1 .25 

0.01696 

787.9 

0.  1822 

1.15 

1 .85 

13. 

1 .2 

0.01541 

790.3 

0.1745 

1.15 

1 .9 

13.2 

1 .2 

0.01443 

789.6 

0. 1676 

1,15 

1 .95 

13.4 

1 .2 

0.01433 

788.9 

0.1678 

1.15 

2. 

13.6 

1 .2 

0.01503 

788.2 

0.174 
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Table  V.  Symbols  used  for  plotting  data  in  Figs.  3-13 
Reference 

Oberlin  (1980,  1982),  Luomajarvi  (1979) 

Bader  (1956),  Borders  (1974) 

Langley  (1976),  Santry  (1984) 

Knudsen  (1980),  Green  (1956),  Chumeinov  (1979) 

Nara 

UCLA,  Lin  (1973) 

DK.  Fontell  (1979) 

Semrad  (1990),  Chu  (1973) 

Sirotinin  (1984),  Kuldeep  (1985) 


symbol 

X 

O 

□ 

K 

0 

X 


\iy 


O 


Nara:  Ishivari  (1984,  1988),  Sakamoto  (1986),  Shiomi  (1986) 

UCLA:  Teasdale(  1 949),  Sonett ( 1 955  ) , Burkig( 1 957  ) , Nakano(l  963  ) 

DK:  Sorensen  (1973),  Andersen  (1967,  1969,  1981) 
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TABLE  VI.  Best  fit  values  of  and  for  57<Z<73,  obtained 
in  a two  parameter  search  for  protons'^  with  T>0.3  MeV  and  a 
with  T>1.6  MeV  with  V *1.25  and  V *1.4,  and  H.  and  from 
Table  III.  a is  defined  by  Eq.  1 1 2 ) , I by  Eq.  ( 1 Obi . If 
V^sO,  H-  is  indeterminate.  The  experimental  errors  a are 
those  given  by  the  authors  (for  some,  an  average  val§e  is 
given).  Authors:  A:  Andersen  et  al.  (1969);  B:  Bader  et 

al.  (1956  );  Bo:  Borders  (1  974  );  C:  Chu  et  al.  ( 1973  );  G: 
Green  et  al.  (1955);  K:  Knudsen  et  al.  (1980);  L:  Luomajftrvi 
(1979);  La:  Langley  and  Blewer  (1976);  N:  Ishiwari  et  al. 
(1988);  0:  Oberlin  et  al.  (1980,  1982);  P:  Lin  et  al.  (1973, 
1974);  S:  Sirotinin  et  al.  (1984);  U:  Teasdale  (1949),  Sonett 
and  MacKenzie  (1955),  Burkig  and  MacKenzie  (1957)  and  Nakano 
et  al.  (1963). 


z 

set 

a 

e 

H3 

^3 

a % 

X 

l./ev 

a 

57 

La 

K 

2 . 

78 

4.7  - 

2.2 

490 

S 

3. 

46.6 

4.47 

0.9 

462 

58 

Ce 

K 

4 . 

75 

4.8 

3.4 

519 

59 

Pr 

K 

3. 

27 

3.75 

2.4 

456 

62 

Sm 

S 

3. 

18 

2.2 

3.4 

548 

64 

Gd 

K 

2. 

48 

3.  1 

1 

582 

S 

3. 

22.5 

2.7 

3.2 

542 

A 

0.6 

0 

0.35 

587 

a 

0 

3. 

79 

4.9 

0.3 

591 

66 

Dy 

K 

2. 

21 

3.3 

0.7 

531 

CL 

P 

4. 

1 9 

3.6 

2.5 

548 

67 

Ho 

K 

2.5 

18.9 

3.58 

2.5 

548 

68 

Er 

K 

3c 

230 

1 0 

2.8 

679 

0 

3. 

33 

1 .3 

1 .0 

691 

La 

2o 

1 0 

2.5 

1 .3 

604 

a 

0 

2. 

18 

2.1 

0.4 

627 

70 

Yb 

K 

2. 

48. 

2.5 

2. 

709 

S 

3. 

18. 

2.15 

2.2 

678 

72 

Hf 

S 

3. 

1 7 

3.5 

1 .6 

600 

39 


73  Ta 


Ta  B 

3. 

10. 

0.6 

0.3 

795 

L 

2.7 

16. 

1 .3 

0.6 

744 

A 

0.6 

10. 

0.6 

0. 

712 

S 

3. 

15.5 

2.2 

1 .4 

706 

U 

- 

- 

0. 

1 .2 

741 

N 

0.3 

14.5 

0.7 

0. 1 

737 

a P 

4. 

21  .5 

1 .5 

0. 

758 

TABLE  Via.  Best  fit  values  of  V for  some  Z obtained  in  a one 

parameter  search  for  protons  with  T>0.3  MeV,  a with  T^1.6  MeV 

with  V -1.25  and  V -1.4,  H and  from  Table  II  and  H.-50 

for  Z<60,  H3-25  for"^  60<Z<72,  H -13  fbr  Z^73.  The  values  ol  V 

are  shown  in  Fig.  4.  and  T are  given.  In  the  last  two 

cols.,  average  I-values , ^q.  (I^b)  are  given,  for  V^-3.85  for 

Z<60,  V =2.3  for  60<Z<72,  V =(Z-46)/25  for  Z>73.  (T^  is  also 

shown.  "^The  latter  I are  also  shown  in  Fig.  5. 

a 


Z 

set 

V 

3 

a % 

X 

I /eV 
a' 

a % 

X 

^a/ 

57 

La 

K 

3.25 

2. 1 

484 

3 . 1 

473 

S 

4.8 

0.9 

462 

2. 

472 

a 

P 

4 

0. 

460 

0. 1 

463 

a 

K 

6.9 

0. 

433 

2. 1 

487 

58 

Ce 

K 

3.05 

3.4 

519 

3.6 

507 

a 

K 

2.4 

510 

59 

Pr 

K 

4.4 

2.9 

491 

3.1 

502 

a 

K 

7.6 

0 

451 

2.9 

528 

62 

Sm 

S 

2.4 

3.6 

558 

3.6 

561 

64 

Gd 

K 

2.45 

1 .6 

555 

1 c8 

561 

S 

2.85 

3.2 

543 

3.8 

560 

A 

0 

0.35 

587 

0.8 

576 

Q, 

0 

2.25 

0.4 

560 

0.4 

558 

a 

K 

1 .8 

0. 

600 

0.5 

573 

a 

C 

548 

66 

Dy 

K 

2.75 

1 .8 

567 

1 .9 

588 

a 

P 

2.2 

0. 1 

609 

0.1 

604 

a 

K 

2.6 

0. 

605 

0.4 

622 

67 

Ho 

K 

3c  1 

3.3 

597 

3.3 

632 

a 

K 

1 .3 

0. 

719 

1 . 1 

655 

40 


68  Er 


K 

2. 

4.3 

635 

4. 1 

620 

0 

1 .3 

1 .4 

682 

2.1 

650 

La 

2.2 

2. 

674 

2 . 1 

671 

a K 

0 

1 .5 

785 

2.7 

633 

a 0 

1 .8 

0.5 

666 

0.8 

640 

a La 

2.8 

0 

608 

0.5 

635 

» K 

2.35 

0.8 

665 

0.7 

667 

S 

2 

3. 

698 

4.1 

684 

a K 

0. 

2.7 

860 

4.4 

691 

72 

Hf 

S 

3.25 

00 

630 

3.2 

671 

a 

C 

683 

73 

Ta 

B 

1 . 1 

0.5 

751 

0.5 

752 

L 

1 . 6 

0.4 

716 

0.8 

751 

A 

0.8 

0.05 

712 

0.2 

707 

S 

2.3 

2.9 

697 

3.7 

760 

U 

0 

1 .2 

741 

1 . 3 

735 

N 

0.5 

0.  1 

738 

0 . 1 

733 

a 

P 

0. 

0.5 

884 

0.7 

787 

82 

Pb 

B 

2.0 

1 .0 

709 

1 .3 

743 

G 

2.9 

1 .7 

652 

2 

753 

A 

1 .4 

0.1 

781 

0.12 

779 

S 

1 .6 

3 

782 

3.0 

794 

U 

0 

0.8 

832 

1 .2 

816 

N 

1 .5 

0.  1 

776 

0.  1 

776 

a 

Bo 

0.7 

0.9 

850 

1 . 

792 

6 
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TABLE  VII.  Energy  loss  measurements  in  thick  absorbers  with 
absolute  energy  measurements:  Bichsel  and  Tschalflr  (1966). 

Initial  energy  T.  and  final  energy  T,  in  The 

experimental  a^sorofer  thicknesses  are  t *0.^385  g/cm'^  for  Pt, 
t -0.5034  g/cm*^  for  Au.  The  theoretical  absorber  thicknesses 
t^(g/cm^)  were  obtained  from  ranges  calculated  with  Eq.  (8), 
and  A*i(t  - t.  ) / t.  . A larger  I-value  would  increase  t . . 
For  a linear^  regression  fit  to  A(t.  ),  the  statistical 
fluctuation  of  A is  ±0.2%. 


^i 

Tf 

^t 

A(%) 

14.364 

4.104 

0.5379 

0.11 

14.416 

4.223 

0.5378 

0,13 

14.538 

4.477 

0.5380 

0. 1 

14.735 

4.909 

0.5370 

0.28 

14.761 

4.964 

0.5368 

0.31 

15.258 

5.952 

0.5354 

0.58 

15.778 

6.881 

0.5352 

0.62 

16.315 

7.785 

0.5349 

0.68 

16.811 

8.585 

0.5343 

0,78 

17.507 

9.651 

0.5340 

0,84 

18.329 

10.866 

0.5329 

1 .06 

19.198 

12.068 

0.5335 

0.95 

20.021 

13.192 

0.5325 

1.12 

20.909 

14.361 

0.5322 

1,19 

21.961 

15.703 

0,5322 

1.18 

23.045 

17.052 

0.5323 

1,17 

24.129 

18.371 

0.5326 

1,11 

25.243 

19.726 

0.5308 

1 .45 

25.433 

19.919 

0.5337 

0,9 

14.011 

4.354 

0.5045 

-0.22 

14.503 

5.339 

0.5042 

-0.15 

15.004 

6.252 

0.5041 

-0.14 

15.513 

7.145 

0.5031 

0.06 

16.16 

8.19 

0.5029 

0.09 

16.821 

9.204 

0.5026 

0.15 

17.631 

10.394 

0.5022 

0.25 

18.459 

11.553 

0.5020 

0,27 

19.306 

12.719 

0.5005 

0.58 

20.171 

13.845 

0.5009 

0.5 

21  .055 

14.976 

0.5009 

0.5 

21 .957 

16.117 

0.4999 

0.71 

22.569 

16.888 

0.4984 

1 ,0 
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TABLE  VIII.  Comparison  <jf  calculated,  S.  , and  measured 
stopping  power,  S (MeV  cm'^  /g)*  The  kinetic  energy,  T,  of 
the  particles  is^  given  in  MeV,  the  authors’  experimental 
accuracy  is  c , and  r(%)  is  the  relative  difference,  Eq.(ll). 


auth 

part 

Z 

T(MeV) 

S ± <7  ( % ) 

X e 

^t 

r% 

Sakamoto 

P 

Ta 

70 

. 4.579  ±0.7 

4.618 

0.9 

1989 

Pt 

fi 

4.508  ±0.7 

4.511 

0.1 

Au 

II 

4.491  ±0.7 

4.518 

0.6 

Pb 

II 

4.388  ±0.7 

4.469 

1 .9 

Chu 

a 

Ta 

2 

366  ±3 

378 

3 

1973 

Au 

II 

358  ±3 

347 

”3 

Pb 

It 

363  ±3 

348 

-4 

Santry 

He^ 

Au 

1 .0 

377  ±5% 

408 

8 

1 981 

1 .2 

363 

380 

5 

1 .4 

348 

357 

2 . 6 

1 . 6 

332 

338 

2 

1 .8 

315 

322 

2 

2-0 

301 

309 

2.6 

Santry 

He"^ 

Au 

1 .0 

390  ±4% 

454 

1 6 

1 980 

1 .2 

380 

424 

1 2 

1 .4 

368 

399 

8 

1 .6 

355 

379 

7 

1 -8 

344 

362 

5 

2.0 

331 

347 

5 

Santry 

He^ 

Au 

3.08 

285  ±8% 

290 

0.4 

1 979 

5.49 

228  ±6% 

223 

-2.2 

Santry 

a 

Au 

3.18 

285  ±4% 

286 

0.4 

1 984 

4.78 

242 

238 

-1  .7 

5.16 

235 

230 

-2.2 

5.49 

228 

223 

-2.2 

5.8 

221 

217 

-1  .8 

6 

221 

214 

-3 

7.69 

207 

1 89 

-9 

Datz 

P 

Au 

2 

43.9  ±0.8% 

45.6 

3.8 

a 

8 

180 

185 

2.9 

Li 

18 

405 

419.5 

3.6 

B 

50 

1 057 

1 1 46 

8 

C 

72 

1429 

1627 

1 4 

N 

98 

1896 

2257 

1 9 

0 

128 

2259 

2815 

25 

P 

3.5 

31.6  ±0.6% 

33.9 

7.4 

a 

1 4 

130.4 

137.2 

5.2 

43 


Li 

31.5 

292 

31  1 

6.7 

B 

87.5 

810 

861 

6.3 

Anthony 

C 

Au 

28.2 

1300  ±2% 

1 509 

1 6 

1 982 

33.6 

1250  ±3% 

1 383 

1 1 

39.6 

1 180  ±3% 

1270 

8 

Ishiwari 

a 

Ta 

8.05 

199.7  ±0.3 

193.9 

"=2.9 

1971,  1977 

28.2 

93.46  ±1 

93.43 

0 

1 978 

28.27 

93.1  ±0.5 

93.28 

0.2 

Au 

28.21 

90.13  ±1 

90.17 

0 

28.28 

91.18  ±0.5 

90.04 

-1  .2 

Takahashi 

a 

Au 

20.6 

no  ±0.7 

109.4 

"0.5 

1983 

27.63 

92.8  ±0.75 

91  .35 

"1  .6 

C 

61  .79 

961  ±0.7 

989 

2.9 

82.9 

839  ±0.7 

828 

-1.3 

a 

Pb 

21  .65 

106  ±0.9 

105.1 

"0.8 

28.59 

89.5  ±0.9 

88.47 

-1.1 

C 

64.94 

982  ±0.7 

951 

-3.2 

85.78 

814  ±0.7 

802.6 

-1  .4 

Bimbot 

0 

Au 

70.24 

1720  ±3 

1 901 

1 0 

1978 

74.88 

1350  ±6 

1835 

36 

75.52 

1350  ±5 

1826 

35 

92.64 

1560  ±4 

1 625 

4 

96.16 

1450  ±8 

1 589 

1 0 

96.96 

1480  ±9 

1 582 

7 

Gauvin 

0 

Ta 

393 

646  ±1.2 

649 

0.5 

1 987 

417 

620  ±1.7 

623 

0.5 

431 

600  ±1.7 

608 

1 .3 

723 

419  ±1.4 

419 

0 

762 

399  ±1 .8 

404 

1 .2 

773 

400  ±2.5 

400 

0 

1418 

263  ±1 .5 

260 

-1.1 

1453 

257  ±2.3 

255 

CO 

o 

1 

Au 

396 

616  ±1.3 

629 

2.1 

421 

600  ±1.7 

602 

0.3 

433 

590  ±1  .7 

591 

0.2 

713 

412  ±1.5 

414 

0.5 

754 

391  ±1.8 

398 

1 .8 

774 

390  ±2.5 

390 

0 

1 350 

263  ±1.1 

263 

0 

1429 

259  ±1  .5 

253 

-2.3 

Bimbot 

Ar 

Ta 

1018 

3050  ±1 

3170 

4 

1986 

1 063 

2920  ±3 

3076 

5 

1 072 

3080  ±4 

3058 

-0.7 

2300 

1800  ±1 .7 

1 794 

-0.3 
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2342 

1 750 

±2.3 

1771 

1 .2 

Au 

935 

3110 

±1.3 

3268 

5 

1 027 

3050 

±1.6 

3067 

0.6 

1 059 

2880 

±3 

3004 

4 

1 073 

2800 

±4 

2978 

6 

1 532 

2330 

± 1 

2329 

0 

1 653 

2150 

±1.5 

2209 

2.7 

1718 

2100 

±2.4 

2151 

2.4 

2157 

1810 

±2.8 

1834 

1 . 3 

2301 

1 730 

±2.3 

1753 

1 .3 

2356 

1 760 

±2.8 

1725 

-2 

2957 

1 440 

±1  .4 

1 472 

2.2 

3028 

1400 

±2.9 

1 448 

3.4 

Ca 

Au 

2919 

1830 

±1.1 

1834 

0.2 

301  0 

1750 

±2.3 

1 796 

2.6 

Schwab 

Ar 

Pt 

2912 

1 366 

±0.6 

1 381 

1 . 1 

1 989 

Au 

2402 

1 556 

±0.  1 

1 581 

1 . 6 

2650 

1451 

±0.3 

1 477 

1 .8 

3005 

1 367 

±0.8 

1 354 

-1 
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TABLE  IX.  Measurements  of  energy  loss  relative  to  A1  or  Cu. 
The  theoretical  values  of  S for  A1  and  Cu  are  given.  The 
Barkas  and  von  Friesen  (1961)  data  for  A1  alloy  have  been 
converted  to  pure  Al^  The  energy  T is  the  average  between  T. 
and  T..  S in  MeV  cm/q,  r given  by  Eq.  (11).  Theoretical  sf 
for  Ai  and  Cu  are  given  for  reference.  Barkas  data  are  for 
pure  Al  ( ICRU  1 984 ) . 


First  author 

Z 

T/MeV 

S % 

X e 

St 

r% 

Bakker 

w 

300 

1 .874 

±1 

1 .858 

-0.9 

Pb 

«t 

1.819 

± 1 

1 .827 

0.4 

U 

H 

1 .736 

±1 

1 .758 

1 .3 

ratio 

Al/Cu 

II 

1.143 

± 1 

1.155 

1 . 1 

Barkas 

Pb 

380 

1.606 

1 .627 

1 .3 

530 

1.391 

1.419 

2 

680 

1 .277 

1 .307 

2.3 

U 

380 

1 . 548 

1 .567 

1 .2 

530 

1 .353 

1 .365 

0.9 

680 

1 .24 

1 .258 

1 . 4 

Al 

380 

2.415 

2.451 

1 . 5 

530 

2.089 

2.116 

1 .3 

680 

1.913 

1 .936 

1 . 2 

Cu 

380 

2.126 

530 

1 .837 

680 

1.681 

Vasilevskii 

Pb 

173.9 

2.496 

±0.8 

2.496 

0 

1 967 

188.5 

2.347 

±0.8 

2.376 

1 .2 

615 

1 .345 

±0.8 

1 .348 

0.2 

650 

1 .332 

±1  .3 

1 .325 

-0.5 

Cu 

173.9 

3.332 

±0.5 

3.32 

-0.3 

188.5 

3.142 

±0.5 

3.154 

0.4 

615 

1 .765 

±0.6 

1 .739 

-1  .5 

650 

1 .734 

±1.0 

1 .706 

-1  .6 

Vasilevskii 

Pb 

1 1 0 

3.314 

±0.8 

3.364 

1 .5 

1 969 

300 

1 .804 

±0.6 

1 .827 

1 .3 

420 

1.531 

±0.6 

1 .557 

1 .7 

600 

1 .343 

±0.5 

1 .359 

1 .2 

Cu 

1 1 0 

4.51 

±0.8 

4.53 

0.4 

300 

2.405 

±0.6 

2.399 

-0.2 

420 

2.029 

±0.6 

2.028 

0 

600 

1 .761 

±0.5 

1 .754 

-0.4 
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TABLE  X.  Comparison  of  experimental  ranges  R to  theoretical 
ranges  R.  (g/cm^).  The  kinetic  energy  the  incident 
particles^  is  T ^(MeV),  Ar  = R - R^ . Except  for  the 
Bloembergen  and  van  Heerden  (19?1)  dara,  the  experimental 
ranges  have  been  corrected  for  multiple  scattering. 


1 author 

Z 

T(MeV) 

R ±<7  % 

X e 

Ar/r  % 

Bichsel 

Au 

9.698 

0.3316 

0.3324 

-0.25 

1 957 

17.549 

0.8685 

0.8667 

0.2 

Bloembergen 

Pb 

62.5 

7.41 

7.375 

0.5 

1 951 

72.1 

9.4 

9.404 

0 

78.9 

11.01 

10.96 

0.4 

82.4 

11.82 

11.8 

0.  1 

89.1 

13.5 

13.48 

0.  1 

- 

95.3 

15.15 

15.12 

0.2 

104.1 

17.58 

17.57 

0.  1 

106.6 

18.28 

18.29 

0 

114.1 

20.7 

20.52 

0.9 

Mather 

Pb 

338.5 

123.25 

123.88 

-0  . 5 

1 951 

339.7 

124.61 

124.58 

0 

Cu 

339.7 

92.69 

93.65 

-1  . 0 

A1 

339.7 

79.42 

80.6 

-1  .5 

Vasilevskii 

Pb 

620 

318.2 

313.4 

1 .5 

1 969 

A1 

620 

- 

206.1 

- 

Barkas 

Pb 

750 

417.7 

413.0 

1 . 1 

1 961 

U 

750 

434.7 

429.2 

1 . 3 

A1 

750 

274.7 

274. 1 

0.2 

Cu 

750 

_ 

316.5 

- 
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TABLE  XI.  Contributions  to  the  stopping  number,  L,  Eqs . (2) 

for  protons  of  energy  T(MeV)  passing  through  a gold  absorber 
with  1 = 790  eV,  i2n  I = 6.672.  L = f(/?)  - Qn  I.  The  shell 
corrections,  c.  and  c are  defiiTed  in  Eq.  (9),  L in  Eq. 
(6),  and  in  Eq.  (/).  The  density  effect  is  included  in 
f{0),  it  amounts  to  0.062  at  100  MeV.  Note  that  a fractional 
change  of  Y%  in  L (or  S)  causes  a change  of  (YL)%  in  the  I- 
value . 


T 

L 

f (/?) 

^i 

-^2 

0.30 

0.5374 

-0.1901 

6.4819 

-0.6270 

0.0398 

0.2338 

0.0935 

0.40 

0.6374 

0.0974 

6.7695 

-0.5398 

0. 1252 

0. 1 968 

0.0713 

0.50 

0.7273 

0.3204 

6.9924 

-0.4724 

0.1800 

0.1721 

0.0577 

0.60 

0.8097 

0.5026 

7.1746 

-0.4175 

0.2163 

0.1543 

0.0484 

0.70 

0.8862 

0.6565 

7.3286 

-0.3714 

0.2408 

0.1407 

0.0417 

0.80 

0.9575 

0.7899 

7.4619 

-0.3317 

0.2573 

0.1299 

0.0366 

0.90 

1 . 0245 

0.9075 

7.5796 

-0.2969 

0.2683 

0.1210 

0.0326 

1 . 00 

1 .0876 

1.0127 

7.6848 

-0.2662 

0.2755 

0.1136 

0.0295 

1.10 

1 .1472 

1 . 1 079 

7.7799 

-0.2386 

0.2798 

0.1073 

0.0268 

1 .20 

1 .2037 

1.1947 

7.8667 

-0.2138 

0.2820 

0.1019 

0.0246 

1 .30 

1 .2574 

1 .2746 

7.9466 

-0.1912 

0.2828 

0.0971 

0.0228 

1 .40 

1 . 3083 

1 .3485 

8.0206 

-0.1703 

0.2823 

0.0929 

0.0212 

1 .50 

1 .3573 

1.4174 

8.0894 

-0.1515 

0.2809 

0.0891 

0.0198 

1 .60 

1 .4041 

1.4817 

8.1538 

-0.1341 

0.2790 

0.0857 

0.0186 

1 .70 

1 .4489 

1 . 5422 

8.2142 

-0.1180 

0.2765 

0.0827 

0.0175 

1 .80 

1 .4920 

1 .5992 

8.2712 

-0. 1 031 

0.2736 

0.0799 

0.0165 

1 .90 

1 . 5334 

1.6531 

8.3251 

-0.0891 

0.2705 

0.0774 

0.0157 

2.00 

1 .5732 

1 .7042 

8.3763 

-0.0761 

0.2672 

0.0750 

0.0149 

2.50 

1 .7533 

1 .9265 

8.5986 

-0.0221 

0.2491 

0.0657 

0.0120 

3.00 

1 .9080 

2.1 080 

8.7801 

0.0181 

0.2309 

0.0589 

0.0100 

3.50 

2.0437 

2.2614 

8.9334 

0.0490 

0.2137 

0.0537 

0.0086 

4.00 

2.1648 

2.3941 

9.0661 

0.0733 

0.1981 

0.0496 

0.0075 

5.00 

2.3737 

2.6156 

9.2876 

0.1081 

0.1712 

0.0434 

0.0060 

10.00 

3.0631 

3.3007 

9.9727 

0. 1 677 

0.0956 

0.0288 

0.0030 

30.00 

4.2076 

4.3686 

1 1 . 0407 

0. 1441 

0.0310 

0.0152 

0.0010 

100.00 

5.4012 

5.4819 

12.1539 

0.0787 

0.0095 

0.0078 

0.0003 
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Table  XII.  Stopping  power  table  for  protons  in  several  elements 
as  a function  of  kinetic  energy  T.  Below  1 MeV,  the  uncertainty 
is  several  percent;  above  1 MeV  it  is  mainly  determined  by  the 
error  in  the  I-value.  If  a linear  interpolation  is  made  for  In  S 
and  In  T,  the  maximum  error  of  interpolated  values  is  0.1%. 


lement : 

La 

Sm 

Er 

Ta 

Au 

Pb 

U 

ensity : 

6.189 

7.49 

9.15 

16.6 

19.32 

11.36 

19.07 

-value : 

474 

561 

650 

734 

790 

779 

841 

T(MeV) 

2 

stopping  power  (MeV  cm  /g) 

La 

Sm 

Er 

Ta 

Au 

Pb 

U 

0.3 

153.79 

125.71 

106.25 

117.29 

103.41 

103.79 

88.36 

0.4 

136.93 

113.95 

98.23 

101.95 

92.02 

92.59 

82.00 

0.5 

123.67 

105.23 

92.06 

91  .86 

84.03 

84.52 

76.32 

0.6 

113.16 

98.12 

86.82 

84.52 

77.99 

78.35 

71,54 

0.7 

104.68 

92.15 

82.24 

78.82 

73.18 

73.43 

67.51 

0.8 

97.707 

87.009 

78.179 

74.208 

69.203 

69.360 

64.055 

0.9 

91 .875 

82.533 

74.553 

70.333 

65.830 

65 . 905 

61.071 

1 . 0 

86.986 

78.624 

71.307 

67.016 

62.910 

62.920 

58.448 

1 .2 

78.882 

72.045 

65.735 

61  .534 

58 . 046 

57.961 

54.032 

1 .4 

72.543 

66.719 

61.136 

57.174 

54.097 

53.959 

50.433 

1 .6 

67.401 

62.304 

57.262 

53.563 

50.820 

50.628 

47,410 

1 .8 

63.113 

58.567 

53.949 

50.499 

48.019 

47.798 

44.825 

2.0 

59.470 

55.344 

51.083 

47.861 

45.586 

45.347 

42.571 

2.5 

52.303 

48.957 

45.324 

42.574 

40.677 

40.409 

38.025 

3.0 

46.964 

44.123 

40.965 

38.557 

36.920 

36.641 

34.528 

3.5 

42.793 

40.317 

37.509 

35.374 

33.926 

33.644 

31  .736 

4.0 

39.419 

37.222 

34.688 

32.770 

31  .470 

31.191 

29.443 

4.5 

36.620 

34.643 

32.334 

30.590 

29.409 

29.135 

27.519 

5.0 

34.254 

32.456 

30.332 

28.734 

27.650 

27 . 382 

25.875 

6.0 

30.451 

28.927 

27.094 

25.720 

24.791 

24.538 

23.206 

7.0 

27.512 

26.191 

24.575 

23.371 

22.553 

22.316 

21.119 

8.0 

25.160 

23.994 

22.549 

21 .475 

20.744 

20.521 

19.435 

9.0 

23.229 

22.185 

20.877 

19.906 

19.249 

19.038 

18.040 

10.0 

21 .609 

20.666 

19.470 

18.585 

17.985 

17.787 

16.863 

12.0 

19.038 

18o247 

17.224 

16.470 

15.960 

15.781 

14.978 

14.0 

17.070 

16.393 

15.503 

14.845 

14.400 

14o239 

13.527 

16.0 

15.527 

14.931 

14.135 

13.551 

13.157 

13.009 

12.370 

18.0 

14.266 

13.736 

13.019 

12.494 

12.140 

12.003 

1 1 .422 

20.0 

13.218 

12.742 

12.089 

11.611 

1 1 .289 

11.162 

10.629 

25.0 

1 1 .232 

10.851 

10.317 

9.926 

9.664 

9.555 

9.112 

30.0 

9.823 

9.506 

9.052 

8.719 

8.498 

8.403 

8.023 

35.0 

8.767 

8.495 

8.099 

7.811 

7.618 

7.533 

7.199 

40.0 

7.943 

7.705 

7.354 

7.098 

6.927 

6.851 

6.552 

50.0 

6.737 

6.546 

6.258 

6.048 

5.909 

5.843 

5.595 

60.0 

5.893 

5.732 

5.487 

5.309 

5.190 

5.133 

4.919 

70.0 

5.267 

5.129 

4.913 

4.757 

4.654 

4.603 

4.415 

80.0 

4.783 

4.661 

4.469 

4.330 

4.237 

4.191 

4.022 

90.0 

4.397 

4.288 

4.114 

3.988 

3.904 

3.862 

3.707 

100.0 

4.082 

3.983 

3.823 

3.708 

3.631 

3.592 

3.450 

eV 
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Table  XIII.  Stopping  power  table  for  alphas  in  several  elements 
as  a function  of  kinetic  energy  T.  Below  2 MeV,  the  uncertainty 
is  several  percent;  above  1 MeV  it  is  mainly  determined  by  the 
error  in  the  I-value.  If  a linear  interpolation  is  made  for  In  S 
and  In  T,  the  maximum  error  of  interpolated  values  is  0.1%. 


Element 

: La 

Sm 

Er 

Ta 

Au 

Pb 

U 

A: 

138.91 

150.35 

167.26 

180.95 

196.97  207.19 

238.03 

I-value 

: 474 

561 

650 

734 

790 

779 

841 

T(MeV) 

stopping  ] 

2 

power  (MeV  cm  /g) 

La 

Sm 

Er 

Ta 

Au 

Pb 

U 

1 .6 

557.42 

466. 1 1 

403.39 

417.59 

378.14 

380.28 

337.98 

1 .8 

529.22 

447.58 

390.41 

395.65 

360.96 

362.94 

326.14 

2.0 

504.1 7 

431.01 

378.50 

377. 1 9 

346.04 

347.86 

315.03 

2.2 

481 .71 

415.85 

367.34 

361 .23 

332.94 

334.50 

304.78 

2.4 

461  .56 

401 .92 

356.82 

3‘47.31 

321.29 

322.61 

295.28 

2.6 

443.45 

389.14 

347.00 

334.94 

310.79 

311.91 

286.51 

2.8 

426.94 

377.24 

337.71 

323.83 

301  .32 

302. 1 8 

278.41 

3.0 

412.05 

366.17 

328.92 

313.78 

292.66 

293.32 

270.88 

3.5 

379.91 

341 .71 

309. 1 1 

292-28 

273.87 

274. 1 2 

254.28 

4.0 

354.18 

321 .03 

291 .82 

274.59 

258.22 

258 . 1 5 

240 . 1 8 

4.5 

332.27 

303.20 

276.65 

259.61 

244.87 

244.54 

228 . 04 

5.0 

313.65 

287.68 

263.31 

246.71 

233.25 

232.76 

217.42 

5.5 

297.56 

274.03 

251.42 

235.45 

223.03 

222.38 

208.05 

6.0 

283.46 

261.91 

240.77 

225.47 

213.92 

213.14 

199.66 

7.0 

259.86 

241.32 

222.50 

208.47 

198.35 

197.43 

185.27 

8.0 

240.73 

224.35 

207.34 

194.48 

185.41 

184.40 

173.26 

10.0 

211.29 

198.01 

183.51 

172.54 

164.98 

163.87 

154.31 

12.0 

189.43 

178.15 

165.55 

155.95 

149.42 

148.28 

139.81 

14.0 

172.39 

162.56 

151.36 

142.85 

137.08 

135.93 

128.29 

16.0 

158.64 

149.92 

139.81 

132.17 

126.99 

125.85 

118.85 

18.0 

147.26 

139.40 

130.20 

123.25 

118.55 

117.43 

110.96 

20.0 

137.65 

130.50 

122.04 

115.67 

111.35 

110.27 

104.24 

25.0 

118.98 

113.15 

106.08 

100.80 

97.22 

96.22 

91  . 04 

30.0 

105.34 

100.42 

94.34 

89.82 

86.74 

85.82 

81.27 

35.0 

94.861 

90.606 

85.267 

81  .309 

78.627 

77.760 

73.698 

40.0 

86.519 

82.775 

78.010 

74.487 

72.101 

71 .304 

67.616 

50.0 

74.003 

70.986 

67.059 

64.164 

62.208 

61.510 

58.407 

60.0 

65.007 

62.456 

59.133 

56.670 

55.009 

54.388 

51 .705 

70.0 

58.181 

56.018 

53.088 

50.947 

49.501 

48.944 

46.574 

80.0 

52.806 

50.914 

48.316 

46.412 

45.135 

44.625 

42.501 

90.0 

48.453 

46.772 

44.424 

42.722 

41  .575 

41.106 

39.180 

100.0 

44.849 

43.336 

41.211 

39.656 

38.611 

38.176 

36.411 
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Figure  captions 


FIG.  1 . Shell  corrections,  calculated  for  gold  with  the 

nonrelativistic  hydrogenic  approximation  (K-shell:  Walske, 

1952;  L-shells:  Bichsel,  1987;  M-shells:  Bichsel,  1983). 

The  scaled  functions  „C..  are  plotted  vs.  a common,  scaled 


abscissa  rj , The  scaling  was 
1.0  at  the  maximum  for  each 
value  at  the  same  value  of 


chosen  to  give  the  ordinate  value 
function  and  to  have  the  maximum 
the  abscissa,  viz.  r/sO-IS.  The 


actual  function,  C^,  can  be  obtained  by  using  rj 


1/ 


and 


^)  =»  where  f^  and 

V 

^1/  are 

K 

given  below 

n 

K 

Is  - - - - 

0.876 

6.03 

1 .04 

2 

LII 

2p  

0.182 

0.999 

1 .472 

2 

MI 

3s 

0.081 

0.616 

0.561 

2 

Mil 

3p  * 

0.072 

0.444 

1 . 04 

2 

MIV 

3d  

0.039 

0.403 

3.524 

4 

The 

functions  for  L^, 

L M 

'*’*TTT 

and  M.. 

are  not  shown  (L^ 

lies 

very 

close  to  M-  ) . 

Note  that  g 

is  not  proportional  to 

the 

number  of  electrons 

n in  the 

subshellls,  and 

only 

approximately  proportional  to  W 


1/ 


Eq.  (4b) 


FIG.  2.  Total  shell  corrections  C/Z  for  protons  with  kinetic 
energy  T in  gold  metal.  Solid  line:  present  values;  dashed- 
dotted  line:  calculation  by  Bonderup  (1967),  using  the  Lenz- 
Jensen  model  of  the  atom;  dotted  line:  experimental  function 
used  by  Andersen  and  Ziegler  (1977);  dashed  line:  Janni 
( 1 982).  Note  that  in  Eq.  (2a)  an  increased  shell  correction 
can  give  approximately  the  same  if  I is  reduced. 


FIG.  3.  Values  of  the  vertical  scaling  factor  for  the 

outermost  shells,  for  best  fits  to  different  experimental  data 

sets  for  protons,  as  a function  of  atomic  number  Z for  Z>73. 

The  parameters  V^=1.25,  ^3“^^  were  used  for  all  Z. 

Resulting  I-values  differ  for  each  set  (Fig.  5).  The 

function  V2=(Z-46)/25  is  shown  as  a solid  line.  It  is  not 

possible  to  assign  uncertainties  derived  from  c to  individual 

e 

values  because  as  well  as  I would  change  with  changes  in 
S^.  Symbols  are  related  to  authors  in  Table  V. 
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FIG.  4.  Same  as  Fig.  3 for  Z<72,  but  data  for  ol  with  T>2  MeV 
(Er,  Z=68,  Oberlin  and  Langley)  were  also  included.  The 
following  parameters  were  used:  V^=1.25,  for  all  Z; 

H^=50  for  Z<60,  H2=25  for  Z>60  (Table  VI).  Resulting  I-values 
differ  for  each  set  of  experimental  data.  No  systematic 
dependence  of  on  Z is  evident,  thus  V^s3.85  was  chosen  as 
an  approximate  average  for  Z<60,  V^s2o3  for  Z>60.  For 
symbols,  see  Table  V. 

FIG.  5.  Values  of  the  Bloch  parameter  b^=I^/Z  for  different 
experimental  data  sets  for  protons  and  a as  a function  of 
atomic  number  Z,  with  V^=1.25,  V2*1.4  for  all  Z,  the  other 
parameters  given  as  a function  of  Z:  for  Z<60,  H^=50,  V^=3.85; 
for  60<Z<72:  H^  = 25,  V^=s2.3,  and  for  Z>72:  V^=(Z- 

46)/25.  The  different  values  b^  at  each  Z express  the 
systematic  differences  in  the  experimental  data  (Table  Via). 
An  uncertainty  of  ±1%  in  gives  an  uncertainty  of  about  ±1% 
in  b^  at  1 MeV,  about  ±3%  at  10  MeV  (Table  X ) . The  unweighted 
mean  value  for  p=18  measurements  in  gold  is  lss(  788±12)  eV,  for 
Pb,  p=8,  I=^(  779±25)  eV. 

FIG.  6.  Comparison  of  experimental  and  theoretical  values  of 

the  proton  and  deuteron  stopping  power  for  gold.  In  order  to 

show  differences  clearly,  the  relative  difference,  r(T)  (Eq. 

(11))  is  given  as  a function  of  the  kinetic  energy  per 

nucleon,  T/M  of  the  particles.  The  theoretical  values  were 

calculated  with  V^=1.25,  V2=1.4,  H^=13,  V^^ ( Z-46 ) /25- 1 . 32 , and 

I =789.9  eV.  Some  data  are  not  shown  here,  but  are  included 
e 

in  Fig.  6a.  Symbols  are  given  in  Table  V;  continuous  lines 
are  used  for  smoothed  data:  dashed-double  dotted:  Bader  at 
al  . (1956);  double-dotted:  Green  et  al.  (1955);  dashed: 
Luomajarvi  (1979);  dashed-dotted:  S?6rensen  and  Andersen 

(1973);  solid  line  above  0.8  MeV:  Andersen  and  Nielsen 

(1981).  The  OL  and  Li  data  by  Andersen  et  al.  ( 1 977  ) between 
1 and  4 MeV/M  are  shown  by  dotted  lines.  For  the  data  by 
Semrad  (1990)  below  0.8  MeV,  a smooth  solid  line  is  shown.  It 
differs  on  the  average  by  less  than  ±0.7%  from  S^,  with  a 
maximum  difference  of  ±1.2%  (see  Fig.  6a).  The  data  by  Green 
et  al.  and  by  Bader  et  al.  were  not  included  in  the  data 
adjustment  of  Section  III.  Thus  their  deviation  is  relatively 
large.  Experimental  uncertainties,  <7^,  given  by  the  authors 
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are  shown  for  only  a few  values.  Negative  values  of  r imply 

values  of  I less  than  I . 

X e 

FIG.  6a.  Comparison  of  experimental  and  theoretical  values 
of  the  stopping  power  for  gold  for  energies  per  nucleon,  T/M, 
between  0.3  and  3 MeV.  Symbols  are  given  in  Table  V. 
Individual  values  of  the  data  by  Semrad  (1990)  are  shown.  a 
data  by  Fontell  and  Luomajflrvi  ( 1 979)  are  shown  as  the  dash- 
dotted  line,  those  by  Matteson  et  al.  (1978)  as  a solid  line 
between  0.35  and  0.55  MeV/u,  for  other  data  see  Fig.  6.  a 
data  by  Lin  et  al.  (1974)  differ  by  less  than  1%  from  those 
of  Fontell  and  LuomajSrvi  and  are  not  shown. 

For  the  Luomojflrvi  and  the  Andersen-Nielsen  data,  similar 
systematic  deviations  of  about  ±0.5%  of  r(T)  are  seen  between 
0.8  and  1.5  MeV.  There  may  be  a problem  in  the  theory  or  in 
the  experiments  or  both  in  this  region.  The  increase  in  r(T) 
for  CL  below  T/M=0.5  MeV  may  be  due  to  a reduced  charge  z of 
the  particles. 

FIG.  7.  Comparison  of  experimental  and  theoretical  values  of 
the  stopping  power  for  tantalum.  The  relative  difference, 
r(T)  in  % (Eq.  (11)),  is  given  as  a function  of  particle 
energy  per  nucleon,  T/M.  The  theoretical  values  were 
calculated  with  H2=13,  V^=1.08,  and  1^-734  eV.  Symbols  are 

given  in  Table  V.  Experimental  error  bars  are  given  at  only  a 
few  values.  Continuous  lines  show  data  smoothed  by  the 

authors.  The  proton  data  by  Bader  et  al.  (1956)  below  0.6 
MeV  are  shown  as  the  dotted  line;  those  of  S^z^rensen  and 
Andersen  (1973)  above  2.25  MeV  as  a solid  line.  The  dashed 
line  shows  the  proton  data  by  Luomajflrvi  (1979).  a-data  by 

Lin  et  al.  ( 1 973  ) are  shown  by  the  solid  line  ending  at 

8.5%.  For  2 MeV  a,  r=(2.9±3)%  for  Chu  et  al.  (1973). 

The  difference  in  r between  the  Luomojarvi  and  the  Andersen 
data  at  neighbouring  energies  (1.5  and  2.25  MeV)  is  about  3%, 
thus  equal  to  the  sum  of  the  experimental  errors  o I 

consider  it  unlikely  that- this  is  a problem  in  the  theory. 

FIG.  8.  Comparison  of  experimental  and  theoretical  values  of 
the  stopping  power  for  lead.  The  relative  difference,  r(T), 
(Eq.  (11))  is  given  as  a function  of  T/M.  The  theoretical 
values  were  calculated  with  H^=13,  V^=1.44,  and  1^=779  eV. 
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For  individual  data  points,  see  symbols  in  Table  V.  Values  of 
r for  the  data  by  Bader  et  al.  (1956),  dashed  line,  and  Green 
et  al.  (1955),  dotted  line,  show  similar  deviation  as  those 
for  gold.  The  data  by  Sorensen  and  Andersen  ( 1 973  ) are  shown 
by  the  solid  line.  The  a-data  by  Borders  (1974)  for 
1 . 2<T(MeV)<1 . 8 are  shown  by  the  solid  line. 

FIGc  9.  Comparison  of  experimental  and  theoretical  stopping 
power  values  for  La,  Ce  and  Pr.  was  calculated  with  the 
parameters  H^^SO,  7^*3. 85.  The  I-values  ( eV ) (Table 
I)  are  shown  next  to  the  chemical  symbol.  Symbols  are  shown 
in  Table  V.  The  solid  lines  represent  the  a-data  by  Knudsen  et 
als  ( 1 980  ),  the  dotted  line  for  La  those  by  Lin  et  al. 
(1973). 

FIG.  10.  Comparison  of  experimental  and  theoretical  stopping 
power  values  for  Sm,  Gd  and  Dy.  is  calculated  with  the 
parameters  H2®25,  and  I^-b  ‘Z  of  Table  I.  Symbols  are 
defined  in  Table  V.  Note  that  Sirolinin  et  al.  (1984)  for  Sm 
and  Gd  gave  two  values  of  S^  at  0.3,  0.4,  0.8  and  1 MeV.  Both 
values  are  shown.  For  Gd,  the  dashed  line  represents  the 
proton  data  given  by  Andersen  et  al.  (1967),  the  dotted  line 
the  OL  data  by  Oberlin  et  al.  ( 1 980  ).  The  solid  lines  show 
the  a-data  by  Knudsen  et  al.  ( 1 980  ),  the  dotted  line  for  Dy 
those  by  Lin  et  al.  (1973). 

FIG.  1 1 . Comparison  of  experimental  and  theoretical  stopping 
power  values  for  Ho,  Er  and  Yb.  S^  is  calculated  with  the 
parameters  and  I^^b^ -Z  of  Table  I.  Symbols  are 
defined  in  Table  V.  The  solid  lines  show  the  a-data  by  Knudsen 
et  al.  ( 1 980  ).  For  Er,  data  by  Oberlin  et  al.  ( 1 982  ) are 
shown  as  smoothed  lines  which  were  obtained  by  making  a three 
parameter  fit  (H^^V^,  I)  to  their  experimental  data;  pi  dotted 
line,  oil  dashed  line.  The  Langley  and  Blewer  ( 1 976  ) data  for 
Er  are  shown  by  the  squares:  empty  for  protons,  full  for  a. 

FIG.  12.  Comparison  of  experimental  and  theoretical  stopping 
power  values  for  Hf  (H^sZS,  v^sa.a),  W and  Ir  (H^^ll,  V2s(z- 
46  )/25).  S^  is  calculated  with  I^  of  Table  I.  For  W,  the 
dotted  line  represents  the  a-data  by  Lin  et  al.  ( 1 973  ),  the 
dashed-dotted  line  those  by  Borders  (1974),  the  dashed  line 
the  proton  data  by  LuomajArvi  (1979),  and  the  solid  line  those 
by  Chumanov  et  al.  ( 1 979  ).  UCLA  data  are  given  for  W and  Ir 
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at  20  and  30  MeV.  For  2 MeV  o. , the  data  by  Chu  et  al.(l  973  ) 
are  plotted  as  a star  at  0.5  MeV. 

FIG.  13.  Comparison  of  experimental  and  theoretical  stopping 
power  values  for  Pt,  Bi  and  U.  is  calculated  with  the 
parameters  H-a13,  V j 2^^  and  of  Table  I.  References 
to  the  symbols  are  given  in  Table  V.  The  proton  data  for  U by 
S^^rensen  and  Andersen  ( 1 973  ) are  shown  by  the  solid  line 
between  2.25  and  18  MeV,  as  are  those  for  Pt  by  Andersen  et 
al.  ( 1 967).  For  protons  in  Bi,  the  data  by  Green  et  al. 
( 1 955  ) are  shown  by  the  solid  line,  those  by  Knudsen  et  al. 
(1980)  by  the  dashed-double  dotted  line,  the  a-data  by  Kuldeep 
and  Jain  (1985)  are  shown  by  the  dotted  line,  those  by  Borders 
( 1 974)  by  the  dashed  line  and  those  by  Knudsen  et  al.  (1980) 
by  the  dasehd-dotted  line.  For  2 MeV  a,  the  value  for  Pt  by 
Chu  et  al.(l973)  is  plotted  as  a star  at  0.5  MeV/u. 

FIG.  14.  Comparison  between  present  calculations  and  some 
other  tabulations.  The  relative  difference  r is  plotted  as  a 
function  of  proton  energy  T,  for  four  elements.  The  solid 
line  represents  the  function  given  by  Janni  (1982),  the  dashed 
line  that  by  Andersen  and  Ziegler  (1977),  the  dotted  line  that 
of  Williamson  et  al  (1966). 


55 


'“’T 


izxt' ' 


/ - 


cittv  i#tn  1 :vn  i n t 


wiO^  ,rt;)4»*  t ^Vj-.  . I‘St^  , ,Tt‘  i9\f  ' ‘ 


fO.'''  10^  &rt^'  .V  mBS0',  fll 

•;|tU 


'Vfcif 


,■/} 


^ W P5 

0^41  «iIodi»rY<* 

U t U /*•  itw6ri«!^.|Nr4t^^«4^^  li4iitfevft4iW4S#'iji4*iP» 

tljpjt  ■■‘-W'-.-.:J:H -I \if,  : Villi'  . ,^ 

» 5.1^1'.  ’©a:?.  < . 

_6  aC^^^fta^^og’jUil»Jlri^i^|^^  "', 

-i##^tiTi»' . <i;.iViit#iWr.»c*ii<4W  baifiioi  ■ ^aijssii^^  -d 

5l,wf:  _ ':Qf  #kie^  rihl^  'niJit ^ 

( p!r€7))  . •' « -lii  . i l*. 'tM  "{.-?  8*?o{l^  bna 

\d  '#;^:t'jv  v#>j  £ ^ 

^.IC^I-  t '■'.  . .i»  "V  -aaSijlsw^mw^lW^^  '' 

.;VaH®|f^r  m,  ^9*?  :^/  '■  ■•:  Oy.^,  ■ C v"'-  ' 

twif  i iHk . ii^ilfcEr 


• '.V  )s^/ 


the  Ly  ;U*U4 

t^o^#  rjy  :;.in  .1,1 


nc* 


1 1 „ 


C-^ji^iur  ;t3^i 


puvfer  ^y{jiix>»»Jt  ^iifr  •'Isilf  Tl| 

^ ' "'  ■’'*  ■m-  "■'  - "■  ■' 


K’^X  fi5  ^ *•  t y ”*"  •^''*  3 f 


- ^*t 

t'V  ■'  C4 

' t>  f I’a  -'^V 


•At  ( :Q:>:^ 

8 iiO.  n.]"  #i 


,V4in^e  -^'h-i^h:  'cJfetH  * rn 


4f4fl  by  « 

l^r  -'“" ‘m  ” '3^1 5»'/..«^i  ‘1  *.n48( , ' , TbV '■  t \ <>yn: 

i'r  -vit'  .'hoiiift/liy , '“’ 


aj 


).ji 


,»i 

ligt ' ' T^^^?J ' j3  ' f <;■  t-  w ^^■:iM»ji 


(f-»  ComJi<ir:.»o6'cJ  «.*p,«ii^itiii!»S" 

ovoi  vaiwoa  ■f.itt  ,«r  ' ili.^*^*ifjpW 

4!5  1 /'iJ  1 .||"  s';,  ■*  - '■■’5' ' - '^'.‘‘'"''.i^'^' 

-aa^hiKl-datt^'j  J I n«,^  .*'^- 

■ .'*»  fii'jjtpn -;^tiife4  wy“  't.ypBja'34¥vljT'^  tpa  isp+w  iia®  ,'-io.*o'  . jj;!-; -s 

> ' ^■.J-’a^!OV''..>^  , asj , (ISTygU  tov  W |jr.:  


•M  - i:t  fct  *'i%'': .,!  f, 

lTn«'  ■>**:o’’'^?  " t't«>73!/  '*•5.0  .^+’  Si 


t^T^b,  tM  yC»kh»c  . Jj 

* tis<l  SoiA*^'  '■■•'ii'tt*  ^«::i..  J 


iJ 


Fig. 


CM 


Bichsel,  5.  Aug.  1990 


C\2 

ob 

« 


z/o 


go 


Bichsel,  5.  Aug.  1990 


zs 


C\2 


o 

zy 


CO 

CO 


CO 

CO 


CO 


C\2 

CO 


O 

CO 


CO 

lO 


CO 

lO 


lO  ^ CQ  C\2  ^ O 


lO 

d 


q 

d 


q q q q 

d d d cd 


(A^rq 


Bichsel,  23.  Aug.  1990 


< 

( 


T/M(MeV) 


Bichsel,  1.  Jan.  1991 


o 

CO 


o 


CO 


CO 

d 


o lo 


o 


iO 


o 


10 

0 

10 

0 

10 

0 

0 


Fig.  9 


Ce,  508  eV  -i 
Z = 58 


1 0.3  13  10 

T/M(MeV) 

Bichsel,  15.  Sep.  1990 


10 

0 

10 

0 

10 

0 

0 


Fig.  10 


1 

c 

— 

) ^ 

I 1 1 I 1 I 1 1 1 1 1 1 1 1 1 

Sm,  561  eV 
r Z = 62 

) o 

u o n 

c 

P c 
) 

- 0 o e 

) -i- 

0.3  1 3 10 


T/M(MeV) 


ir;  100^^ 


10 

0 

10 

0 

10 

0 

0 


Fig.  11 


□ -L 


Er,  650  eV 
Z-68 


B □ 


-E^ 


O 


C) 


$ o 
o 


Yb,  676  eV 
Z-70 


■Q...-Q a 


o 


6 ^ 

-L 

~ 

\ ^ 

V 1 111 

_ , , ,J 



1 III 

1 till 

0.3 

1 

3 

10 

T/M(MeV) 


Bichsel,  15.  Sep.  1990 


Fig.  12 


Bichsel,  15  Sep.  1990 


Fig.  13 


T/M  (MeV) 


Bichsel,  1.  Jan.  1991 


Fig.  14 


Bichsel,  15.  Sept.  1990 


NIST.114A  U.S.  DEPARTMENT  OF  COMMERCE 

(REV.  3-89)  NATIONAL  INSTITUTE  OF  STANDARDS  AND  TECHNOLOGY 

BIBLIOGRAPHIC  DATA  SHEET 

1.  PUBUCATION  OR  REPORT  NUMBER 

NISTIR  4550 

2.  PERFORMING  ORGANIZATION  REPORT  NUMBER 

3.  PUBUCATION  DATE 

APRIL  1991 

4.  TITLE  AND  SUBTITLE 

Stopping  Power  of  Fast  Charged  Particles  in  Heavy  Elements 

5.  AUTHOR(S) 

Hans  Bichsel 

6.  PERFORMING  ORGANIZATION  (IF  JOINT  OR  OTHER  THAN  NIST,  SEE  INSTRUCTIONS) 

U.S.  DEPARTMENT  OF  COMMERCE 

NATIONAL  INSTITUTE  OF  STANDARDS  AND  TECHNOLOGY 

GAITHERSBURG,  MD  20899 

7.  CONTRACT/GRANT  NUMBER 

8.  TYPE  OF  REPORT  AND  PERIOD  COVERED 

9.  SPONSORING  ORGANIZATION  NAME  AND  COMPLETE  ADDRESS  (STREET,  CITY,  STATE,  ZIP) 

10.  SUPPLEMENTARY  NOTES 


DOCUMENT  DESCRIBES  A COMPUTER  PROGRAM;  SF-185,  RPS  SOFTWARE  SUMMARY,  IS  ATTACHED. 

11.  ABSTRACT  (A  200-WORD  OR  LESS  FACTUAL  SUMMARY  OF  MOST  SIGNIFICANT  INFORMATION.  IF  DOCUMENT  INCLUDES  A SIGNIFICANT  BIBUOGRAPHY  OR 
LITERATURE  SURVEY,  MENTION  IT  HERE.) 

The  stopping  power  formula  from  Bethe's  theory  contains  terms  which  are 
known  only  approximately  and  must  be  estimated  with  the  use  of  experimental 
data.  These  terms  include  a material  constant,  the  mean  excitation  energy  of 
the  medium,  and  the  shell-,  Bloch-  and  Barkas - corrections . In  an  analysis  of 
measured  proton  and  alpha-particle  stopping  powers  and  ranges,  modifying 
parameters  have  been  introduced  into  these  corrections , and  the  mean 
excitatibn  energy  was  simultaneously  adjusted,  so  as  to  get  the  closest 
possible  agreement  with  experimental  results.  Such  an  analysis  is  reported 
here  for  elements  with  atomic  numbers  Z > 57.  The  modification  parameters 
introduced  for  the  shell  corrections  have  a simple  relation  to  atomic  energy 
levels.  The  Bethe  theory  with  the  adopted  mean  excitation  energies  and 
proposed  adjustments  of  the  shell-  and  Barkas -corrections  predicts  stopping 
powers  that  are  in  close  agreement  with  experimental  values , within  the 
experimental  uncertainties.  This  agreement  was  obtained  for  protons  with 
kinetic  energies  above  about  0.5  MeV,  and  for  heavier  ions  of  charge  z at 
energies  above  (z  - 1.5)  MeV/u. 


12.  KEY  WORDS  (6  TO  12  ENTRIES;  ALPHABETICAL  ORDER;  CAPITAUZE  ONLY  PROPER  NAMES;  AND  SEPARATE  KEY  WORDS  BY  SEMICOLONS) 


Bethe  theory;  charged  particles;  heavy  elements;  mean  excitation  energies; 
protons;  shell  corrections;  stopping  power 


13.  AVAILABILITY 

14.  NUMBER  OF  PRINTED  PAGES 

X 

UNUMITEO 

73 

FOR  OFFICIAL  DISTRIBUTION.  DO  NOT  RELEASE  TO  NATIONAL  TECHNICAL  INFORMATION  SERVICE  (NTIS). 

ORDER  FROM  SUPERINTENDENT  OF  DOCUMENTS,  U.S.  GOVERNMENT  PRINTING  OFFICE, 

WASHINGTON,  DC  20402. 

IS.  PRICE 

A04 

X 

ORDER  FROM  NATIONAL  TECHNICAL  INFORMATION  SERVICE  (NTIS),  SPRIflGFIELD,  VA  22161. 

ELECTRONIC  FORM 


